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In a wide class of holographic models, like the one proposed by Sakai and Sugimoto, 
baryons can be approximated by instantons of non-abelian gauge fields that live on 
the world-volume of flavor D-branes. In the leading order, those are just the Yang-Mills 
instantons, whose solutions can be constructed from the celebrated ADHM construction. 
This fact can be used to study various properties of baryons in the holographic limit. 
In particular, one can attempt to construct a holographic description of the cold dense 
nuclear matter phase of baryons. It can be argued that holographic baryons in such a 
regime are necessarily in a solid crystalline phase. In this review we summarize the known 
results on the construction and phases of crystals of the holographic baryons. 

Keywords : cold dense nuclear matter; Sakai-Sugimoto model; holographic baryons; in¬ 
stantons. 


1. Introduction 

The 1980’s witnessed a revival of interest in the Skyrme model of baryons 1 At the 
time there was a growing evidence that the correct low-energy description of baryons 
is that of solitons of meson fields!^ The chiral Lagrangian, unknown at the time of 
the Skyrme’s original work, provided a natural framework for the realization of this 
ideaP Starting with the work of Adkins, Nappi and Witten^, many static properties 
of baryons were computed using the Skyrme model. The approach produced many 


* Prepared for the special issue Skyrmion of the IJMPB, ed. by M. Rho and I. Zahed. 
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new exciting results (for reviews see e.g. Refs.fSHJ), though remained far from being 
a systematic accurate method to describe baryons. 

A branch of applications of the Skyrme model, which partially motivated our 
work and the current review, emerged in the analysis of the properties of cold 
dense baryonic matter. In 1985 Klebanov® analyzed skyrmions arranged in a three- 
dimensional simple cubic lattice and found a ground state of such a system. How¬ 
ever, the true ground state of skyrmions at finite density, is not the simple cubic, 
but rather the face-centered cubic lattice. This was later found by Kugler and 
Shtrikman.^ E3 Even before that it was noticed by Goldhaber and Mantorfm, that 
at large density — essentially when the skyrmions start to overlap — skyrmion 
lattices undergo a transition to configurations with a higher degree of symmetry. 
In the new state, the appropriate description of the system seemed to be not the 
one of the skyrmion, but rather of the half-skyrmion lattice. In particular, the ele¬ 
mentary Wigner-Seitz cell of the high density configuration contained a half of the 
skyrmion’s topological charge. Kugler and Shtrikman found that the low-energy 
ground state, the face-centered cubic lattice of skyrmions, at high densities turn 
into a simple cubic lattice of half-skyrmions. 

An important property of the half-skyrmion phase is the partial restoration of 
chiral symmetry, which occurs when skyrmion fields are averaged over an elementary 
cell. Such behavior can be compared with baryons at finite density. The N c = 3 QCD 
at the chemical potential slightly larger than the baryon mass is believed to be in a 
deconfined color superconductor phase with the chiral symmetry restored. For large 
N c there is no color superconductivity and no deconfinement, for reasonably large 
values of the quark chemical potential gL q Aqcd- In this regime, the phase is 
some kind of a quark liquid, with quark like excitations in the interior of the Fermi 
sea, but with hadronic ones close to the Fermi surface!^ The interior of the Fermi 
sea is chirally symmetric, but the symmetry is broken close to the surface by the 
chiral density waves. In other words the chiral symmetry is locked with translational 
symmetry, exactly like in the case of skyrmions. 

Another part of the motivation to study baryonic crystals came from the ad¬ 
vances in the holographic correspondence, or simply holography^in the 2000’s. 
Shortly after the original proposal of the AdS/CFT correspondence,^ it was under¬ 
stood how to put baryons on the holographic setting, e.g. in Refs. H4UI6I However, 
the most popular setup to study holographic hadron physics came a bit later, with a 
proposal of Sakai and Suginroto in 2004!^ Remarkably, the latter model naturally 
incorporates in it the original Skyrme model with some modifications. Specifically, 
in terms of the pion fields, one can present the effective action as the chiral La- 
grangian plus the Skyrme term, plus an infinite series of the couplings to vector 
meson fields. 

For a reader less acquainted with the developments in the holographic corre¬ 
spondence, let us list the main milestones as far as applications to hadron physics 
is concerned. The correspondence provides a very powerful and useful tool to study 
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strongly coupled systems by mapping them to a weakly coupled string theory, equiv¬ 
alently low-curvature gravity. The most natural physical system in which it could 
be and was implemented is QCD and hadron physics. 

It is well known that the gauge theory dual of the type IIB string theory on 
the AdSfj x S 5 space is the AT = 4 super Yang-Mills (SYM) theory, which is con¬ 
formal, with the maximal amount of supersymmetries. In order to study purely 
non-supersymmetric YM theory one has to deform the gauge theory so that all 
supersymmetries are broken and the system turn into a confining rather than con¬ 
formal phase. Associated with such a deformation the background string theory has 
to be deformed accordingly. 

A prototype holographic model of this nature is the so called Witten’s model, 
where one starts with the near horizon limit of the background of large number N c 
of D4 branes (instead of D3 branes in the AdS$ x S 5 scenario) and compactifies 
one space coordinate in such a way that together with the radial direction it has 
a cigar-like shape. Imposing anti-periodic boundary conditions for the fermions 
breaks all the supersymmetries so that in the limit of a small compactification 
radius the geometry is dual to the low-energy regime of the four-dimensional YM 
theory contaminated with Kaluza-Klein modes. 

This background admits a Wilson loop with an area law, namely confining, 
behavior as well as a gapped glueball spectrum. To incorporate quark degrees of 
freedom one introduces additional Nf “flavor” branes, typically in the probe Nf <C 
N c regime. This prescription was originally introduced by Karch and Katz,® and 
first implemented for a system in the Coulomb phase. It was late r co nsidered in 
a confining background in Ref. [2U1 and finally, Sakai and Sugimotcff^ proposed a 
model that up to date provides the most useful holographic dual model of a system 
in the same “universality class” as QCD with Nf << N c flavors. 

The model in based on extending the Witten’s model by placing stacks of Nf 
D8 and anti-D8 flavor branes that merge with each other at the tip of the cigar, 
so that the construction exhibits a “geometrical spontaneous breakdown” of chiral 
symmetry. This model was later generalized by Aharony et a/RU by allowing the 
flavor branes to merge at any point along the radial direction and not necessarily 
at the tip of the cigar. (See a more detailed discussion below.) 

A Wilson line is described in holography by a string that stretches out and 
ends on the boundary of the holographic background. The endpoints of the string 
can be viewed as an infinitely heavy quark and anti-quark pair. In this picture a 
meson takes the form of a string whose endpoints are flavored and of finite mass. 
In holography this maps into a string that starts and ends on probe flavor branes. 
In the limit of large string tension and small string coupling, which corresponds in 
the gauge theory side to the limits of large number of colors N c and large’t Hooft 
coupling A = jy M lV cl one can describe the degrees of freedom of the meson by those 
of fields (scalars or vectors) that reside on the flavor branesP^ On the other hand 
a description of realistic mesons, and in particular their Regge spectrum, requires 
the use of the stringy holographic picture!^ 
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What is the holographic dual of a baryon? Since a quark corresponds to an 
endpoint of a string, a baryon should include a structure to which N c strings are 
connected. It was found out that in the context of the AdS 5 x S 5 string theory such 
a construction takes the form of a “baryonic vertex”, which is a D5 brane that wraps 
the S 5 with N c strings that connect it to the boundaryP3 This in fact corresponds 
to an external baryon with infinitely heavy quarks. Such construc tions o f stringy 
baryons were proposed also in the context of confining backgrounds 2 D. 

In a similar manner to the passage from a Wilson line to a meson, to get a dy¬ 
namical baryon, and not an external one, the strings have to connect the baryonic 
vertex to flavor branes and not to the boundary. It was shown in Ref. [24] that in 
the generalized Sakai-Sugimoto (gSS) model the baryonic vertex is immersed in the 
flavor brane itself. Recall that in the gSS model the flavor branes are D 8 branes 
that wrap an S 4 and hence the baryonic vertex, which in this model is a D4 brane 
wrapping the same cycle, is a soliton to the five dimensional flavor gauge theory 
that resides on the flavor branes. It is also an instanton in the Euclidean four di¬ 
mensions spanned by the three ordinary space directions and the radial holographic 
direction.251 Thus, in the large N c and large A limit, baryons can be described as 
flavor instantons. 

Similarly to the Skyrme model, the Sakai-Sugimoto model can be analyzed to 
extract the static properties of baryons. In certain aspects it may do better then 
the Skyrme model, but overall its predictive power is comparable!^--' What is inter¬ 
esting about the holographic realization, is that static baryons can be described by 
solitons (instantons) of the non-abelian gauge (Yang-Mills) fields. Curiously, simi¬ 
lar approximate description was ingenuously anticipated by Atiyah and Manton for 
skyrmions somewhat a decade before holography. ' 

The main problems one encounters when describing baryons holographically are 
the following. First, the size of the baryon turns out to be of the order 1/y/X and 
hence in the A —> 00 limit it becomes comparable to the string scale and stringy 
corrections are not necessarily negligible. Second, there are two different scales, off 
by a factor ~ 2, one has to use to describe the mesonic and baryonic spectra. The 
latter problem, but not the former, can be overcome when one uses the gSS rather 
than the original SS model.^2 

Another severe problem in using the SS model is the fact that the interaction 
between the baryons of this model is repulsive at any separation distance. For long 
distances between the baryons the repulsion is due to the fact that the lightest 
isoscalar vector, whose exchange yields repulsion, is lighter than the lightest scalar 
that yields attraction. In the near and intermediate zones the interaction between 
two instantons of the SS model is purely repulsive. By using the gSS instead of the 
SS model, the severeness of the problem can be reduced. As was shown in Ref. [283 
in the gSS (and not in the SS) model there is, in addition to the repulsive force, 
also an attractive one due to an interaction of the instantons with a scalar field that 
associates with the fluctuation of the embedding, though the ratio of the attractive 
to repulsive potential can never exceed 1/9. Thus in both the SS and the gSS nuclei 
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will not be formed. 

In Ref. [ 37 ] it was shown that there is another holographic model^ with a dom¬ 
inance of the attraction at long distances, but at the same time, a tiny ~ 1.7% 
binding energy. In that model the lightest scalar is in fact a pseudo-Goldstone bo¬ 
son associated with the spontaneous breaking of the scale symmetry. This meson 
can be made to be lighter, but not much lighter, than the lightest isoscalar vector 
and hence the emerging system is that of small attraction dominance.— The fact 
that in the gSS the interaction between the baryons is for any separation distance 
repulsive will be important, but not crucial for this review, as we shall see in more 
detail later. 

Motivated by the interesting behavior of skyrmions at high density and being 
equipped with the new methodology of holography, in Refs. l29l and l30l the authors 
of this review decided to look at the problem of the cold nuclear matter from the 
holographic perspective. Conversely to the real nu clea r matter, in the large N c limit 
of holography the latter is always a solid (crystal) J2S! Thus, it appeared interesting 
to look at the realization of the crystal of baryons in the instanton description of 
the Sakai Sugimoto-like setup. Moreover, related studies appeared in papers IMU351 
In particular, Refs. [32] and [33] have analyzed the phase diagram of configurations 
with finite baryon density, approximating the density of instantons by a mean value. 
Refs. [34] and [351 searched for the manifestation of the half-skyrmion transition in 
a holographic setting. The research program, to be described here, of studying the 
instanton crystals as holographic duals of the cold nuclear matter stemmed from 
these and related results. 

In brief, the program assumes construction of instanton solutions that may ap¬ 
proximate the crystals of the holographic baryons. We use the Sakai-Sugimoto model 
as a prototypical example, but any other model where the baryons are realized as 
instantons would be just as good for our purposes, for example the model of Ref.foITl 
for some seven-brane geometriesor the AdS .5 x S 1 model of Ref. [38] (the baryons 
of that model are studied in Sec. 6 of Ref. [24}. For this large class of models one can 
introduce a general energy functional that will need to be minimized over various 
instanton configurations to find the true ground state for a given set of external 
parameters. 

Holographic study implies certain limits for the parameters of the underlying 
gauge theory, such as large number of colors N c and large’t Hooft coupling A. On 
one hand this simplifies the problem and allows to use perturbation theory in the 
construction of instanton solutions. One the other hand it introduces new problems, 
such as additional suppression of the baryon radii r ~ which is problematic 

not only for the phenomenology of the model, but also for its general consistency. 
In the leading order (LO) of the perturbation theory the solutions are just the 
well-known instantons of the SU ( N ) Yang-Mills. However all these instantons are 
LO degenerate and one has to minimize the energy in the next-to-leading order 
(NLO), selecting from a large moduli space of the LO solutions. The problem could 
be simplified by appropriately reducing the moduli space, but the experience shows 
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that the our naive intuition does not always work well in selecting the appropriate 
configurations. 

The review is supposed to summarize our findings and expectations in the re¬ 
alization of this program. We will assume that the cold nuclear matter phase of 
the holographic baryons exist and is realized by an appropriate configuration of the 
Yang-Mills instantons. We will discuss examples of instanton configurations and try 
to select an appropriate ones to describe the ground states. 

The review is organized as follows. In section [2] we review how baryons can 
be realized in the Sakai-Sugimoto and a class of similar holographic models. In 
section EH we review the model of Sakai-Sugimoto, while in section f2. 2 1 we discuss 
baryons, their properties and further adaptation and generalization of the model, 
suitable for the subsequent study. 

In section [3] we discuss the consequences of the holographic limit (large number 
of colors and large ’t Hooft coupling) for the physics of hadrons at finite density. 
In subsections IQ and 13.21 we consider the general implications of these limits. In 
particular, we explain why we expect to find baryon crystals at large N c and outline 
the main features of such a phase. In section 13.31 we review what the holographic 
limits mean specifically in the context of the Sakai-Sugimoto and similar models. We 
outline the main steps needed for the instanton description of the baryon crystal. 
In section l3.4l we illustrate our expectations in a simplest example using the lattice 
of point charges. This example summarizes the behavior of the holographic lattices, 
when the density is increased, and demonstrates what we believe is the onset of the 
transition to the quark liquid phase. 

In section 01 we review available exact solutions for the instanton lattices. These 
are only known for the (quasi-) one-dimensional configurations, a straight periodic 
chain and a (abelian) zigzag. Naively the zigzag is the first configuration to which 
the straight chain can break, when the pressure is increased. However, only the 
straight chain can be realized as a ground state in our setup. For the efficient study 
of more complicated configurations it seems necessary to rely on approximations. 
This is implemented in section [5] 

In section lSTl the two-body force approximation is discussed. This approximation 
can be used in our setup and it allows numerical simulations to be used for the 
search of available energy minima. It is employed in section 15.21 where the possible 
quasi-one-dinrensional configurations are considered. It is assumed that at densities 
not too high possible configurations are either linear arrays (straight chains) of 
instantons, or zigzag-like configurations. Section 15.31 summarizes the results of a 
more recent study in progress of the two-dimensional lattices, both infinite and 
finite size. It demonstrates that if the zigzag assumption of Sec. 15.21 is relaxed, there 
exist other configurations that can be favored over the (abelian) zigzag. 

Section [6] summarizes the main results and problems. 
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2. Holographic QCD 
2.1. Sakai-Sugimoto Model 

Some gauge theories, e.g. AT = 4 supersymmetric Yang-Mills (SYM) theory, have 
exact holographic duals, where both sides of the duality follow as IR limits of the 
same string-theoretical construction, while all the undesirable degrees of freedom 
are superheavy. This is unfortunately not the case of Quantum Chromo Dynamics 
(QCD), not even of its large N c limit: it either does not have an exact holographic 
limit, or we have not found it yet. Instead, there is a large number of “holographic 
QCD” (hQCD) models, which are rather dual to some QCD-like theories with a lot 
of extra stuff that the real QCD does not have. Our hope is that such models are 
good enough to get qualitative understanding, though they cannot probably be too 
accurate to make useful numerical predictions. 

One of the most popular hQCD models on the market is the Sakai-Sugimoto 
modeff^, which we are going to review momentarily. The predictions of this model 
are at best qualitative, although it can successfully compete with conventional 
hadron physics models like Skyrrne model as far as certain quantitative predictions 
are concerned. 

The construction of the Sakai-Sugimoto model starts with N c coincident D4 
branes, spanning the Minkowski space times a circle of radius R with antiperiodic 
boundary conditions for the fermions. Such a configuration breaks the H = 4 su¬ 
persymmetry down to J\f = 0*. For weak ’t Hooft coupling A = gy M Yc ~ g s -/V c 1 
(field theory limit, g s — string coupling), the open strings between the branes 
give rise to the gluons of the pure U(N C ) Yang-Mills (YM) theory. However in 
the holographic limit A /?> 1 everything happens right at the Kaluza-Klein scale 
Aqcd ~ Mkk = 1/i?, so the YM glueballs end up with similar O(Mkk) masses to 
a lot of non-YM stuff. This scenario is known as the Witten’s model.^® 

On the gravity side of the duality, at A> 1, the D4 branes merge into a black 
brane which warps the 10D metric. All we see outside the horizon is a warped space- 
time geometry and the Ramond-Ramond flux induced by the conserved charge 
of the D-branes. Specifically we have a warped product of R 3,1 Minkowski space, 
the S 4 sphere (originally surrounding the D4 branes), and a two-dimensional cigar 
spanning the radial direction perpendicular to the branes and the S' 1 circle. The 
radial coordinate, here denoted u, runs from ua > 0 to infinity. At ua the S' 1 shrinks 
to a point, hence the cigar. 

Altogether, we have warped metric, the four-form flux, and the running dilaton 
according to 


ds 2 = 


3/2 


Rda) 


—dt 2 + 


Sijdx l dx 3 


f(u)da 


+ 


Rd4 


3/2 


du 2 

J(u) 


2 dtti 


F 4 = Sni^Nc x vol(S 4 ), 



e 


9s 


( 1 ) 
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where X 4 is the coordinate along the S' 1 or the polar angle on the cigar, 

R 3 D4 = *g,%N c , /( u ) = l-(^) 3 , 


( 2 ) 


l s = vcd is the string length scale. The u\ — the minimal value of the radial 
coordinate u at the tip of the cigar — is related to the original radius R of the S 1 
circle as 


2irR 


4^ ( Rd±\ ! 
3 V UK J 


( 3 ) 


The same radius R also controls the 4D Yang-Mills coupling and hence the ’t Hooft’s 
coupling A. Analytically continuing from A 1 to A 3> 1, we have 

27t g s i, 


^ — 5Id N c — 


95D_ 

2tt R 


N r . = 


R 


■AT 


( 4 ) 


To add the flavor degrees of freedom to the model, Sakai and Sugimoto added 
Nf D8 and Nf anti-D8 (flavor) branes. The D8 (anti-D8) branes span all space 
coordinates except the X 4 coordinate along the S 1 circle. At weak A, the open 
strings connecting the color branes to the flavor branes give rise to the quarks and 
the anti-quarks. The massless quarks are localized at the intersections of D8 and 
D4 branes; likewise, at the intersections of the anti-D8 branes and the D4 branes we 
get massless antiquarks. The open strings between the flavor branes yield 7V| vector 
and scalar fields living on those branes. The four-dimensional modes of these vector 
and scalar fields are dual to the QCD mesons. In the meantime the YM instantons 
of the vector fields are dual to the QCD baryons, see Sec. 12.21 for more details. 

On the holographic side of the duality (N c —> 00 , A —> 00 ) the exact solution for 
the flavor branes interacting with the warped metric and fluxes is not known, but 
for Nf < N c and g s Nf -C 1 we may use the probe approximation: the flavor branes 
seek the lowest-action configuration in the background metric CD, while their back- 
reaction upon the metric is neglected. The flavor branes remains D-branes, rather 
than merge into a black brane of their own. 

Consequently, at low temperatures (below the deconfinement transition]^, the 
flavor branes span a product of the Minkowski space, four-sphere S 4 and a one¬ 
dimensional curve on the cigar; the exact shape of this curve follows from minimizing 
the branes’ action, but its topology follows from the cigar itself: since the D8 and 
the anti-D8 branes cannot continue all the way to the origin of the geometry u = 0, 
they must reconnect to each other and form a U-shaped configuration as shown on 
figure [Q 

The reconnection is the geometric realization of the spontaneous chiral symmetry 
breaking: the separate stacks of Nf D8 and Nf anti-D8 probe branes give rise to the 
U(Nf)L x U(Nf)n gauge symmetry, which corresponds to the U{Nf)i, x U(Nf)n 


a At higher temperatures — above the deconfining transition — the background metric has dif¬ 
ferent topology, and the flavor branes also have different shapes, see Ref. 1211 for details. 
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Fig. 1. The figure on the right is the generalized non-antipodal configuration. The figure on the 
left describes the limiting antipodal case L = nR, where the branes connect at uq = u\ . 


global chiral symmetry in the dual four-dimensional theory. But when the D8 and 
anti-D8 reconnect, only a single stack of Nf U-shaped branes remains and hence 
only one unbroken U(Nf) symmetry. Thus, the chiral symmetry is spontaneously 
broken, 


U(N f ) L x U(N f ) R -»■ U(N f ). (5) 

As shown on figure Q3 the U-shaped profiles of reconnected branes depend on 
an additional parameter — the asymptotic separation L of the D8 and the anti-D8 
branes along the S 1 circle for u —> oo. For L = nR the branes form the antipodal 
configuration in which the branes remain at opposite points on the circle for all 
it: this is the original configuration of Sakai and Sugimoto. In the more general 
version of the modeb^ we allow for the L < irR non-antipodal configurations. In 
such configurations, the distance between the branes in the X 4 direction depends 
on u - it becomes smaller for smaller u — and eventually the branes reconnect at 
uq before they reach the bottom of the cigar. The £ = Uq/u\ ratio may be used to 
parameterize the non-antipodal configurations instead of the L/R. 

The ( or L/R parameter of the Sakai-Sugimoto model does not correspond to 
any known adjustable parameters of the real-life QCD. Unfortunately, this parame¬ 
ter affects many physical properties of the model. For example, for (L/R) > 0.97 the 
deconfinement and the restoration of chiral symmetry happen at the same temper¬ 
ature, but for (L/R) < 0.97 they happen at different temperatures and the model 
has an intermediate deconfined but chirally broken phase'll Also, in the antipodal 
model the central nuclear forces are purely repulsive, while the non-antipodal mod¬ 
els give rise to both repulsive and attractive nuclear forcefPSl, though the net force 
remains repulsive at all distances. 

The low-energy dynamics of the flavor degrees of freedom living on the D8 branes 
is governed by the effective action comprised of the Dirac-Born- Infeld (DBI) and 
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Chern-Simons (CS) terms, 


S = S'dbi + Scs ■ 


( 6 ) 


The DBI action is 

-S'dbi = -Tg J d 9 x e -0 Str det (g mn + J" m „)) > ( 7 ) 

D8+m 

where Tg = {2 tt)~°£~ 9 is the D8-brane tension, g mn is the nine-dimensional induced 
metric on the branes, F mn is the U(Nf) gauge field strength, and Str denotes the 
symmetrized trace over the flavor indices. The lowercase arabic letters m, n denote 
the nine coordinates of the D8 branes’ world-volume. 

In the limit of fixed brane geometry and weak gauge fields, the DBI action 
reduces to that of Yang-Mills, 

S'dbi [-T 7 ] = const + Sym[F\ + 0(T 4 ). (8) 


Furthermore, the low-energy field modes we are interested in are constant along the 
sphere S 4 and only contain the components of vector fields perpendicular to S 4 . 
Thus we are going to dimensionally reduce the flavor gauge theory down to five 
dimensions: the four Minkowski dimensions x 0,1 ’ 2,3 , plus one coordinate z along the 
U-shaped line on the cigar. We find it convenient to choose a particular z coordinate 
that makes the five-dimensional metric conformal 



In the x M = (x°, x 1 , x 2 , x 3 , z) coordinates, the five-dimensional YM action for 
the flavor gauge fields becomes 


Sym 


f d 4 x / dz 


2 9ym ( z ) 


tr 0 ^mn) 


where 

1 _ N c \Mkk u(z) 

2 5ym(- z ) 216 tt 3 X u a 

Near the bottom of the U-shaped flavor branes one finds 


( 10 ) 


( 11 ) 


1 

2 .9ym( 2 ) 


N c XMkk 

2167T 3 



8C 3 - 5 


M^ k z 2 + 0(M 4 K 



( 12 ) 


The Chern-Simons term arises from the coupling of the gauge fields on the D8 
brane to the bulk Ramond-Ramond field. In nine dimensions 


Scs = T s J C 3 Atre w ^, 

D8+D8 


where F 4 = dC 3 . 


(13) 
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After integrating over the S 4 and dimensionally reducing to five dimensions, the CS 
term becomes 

Scs = < 14 > 

In a particularly interesting case of two flavors, it is convenient to separate the U (2) 
gauge fields Am into their SU(2 ) components Am and the C/(l) components Am- 
In terms of these components, the CS action becomes 

+ 9< 15 > 


We shall see in a moment that the baryons and the multi-baryon systems have 
strong self-dual SU (2) magnetic fields 0 Thanks to the first term in this CS 
action, the instanton number density, 


I(x,z) 


1 

32t r 2 


rpa rpa 


(16) 


acts as electric charge density for the abelian field Aq ; the net electric charge of an 
instanton is Q e \ = N c / 2. 

Besides the U(Nf) gauge fields, the effective low-energy five-dimensional theory 
also contains the scalar fields d ,a (x,z) describing the small fluctuations of the D8 
branes in the transverse directions. For Nf branes, the scalars form the adjoint 
multiplet of the U(Nf) gauge symmetry. The action for the scalar fields follows 
from the DBI action for the induced metric g mn of the fluctuating branes. For the 
$(*, 2 ) fields normalized to have similar kinetic energies to the vector fields, the 
scalar action looks like 


Scalar 


+ Jd A x Jdz C(z) x tr($J r A/A tF mn ) + ••• . 


(17) 


The details of the scalar potential V($) = m 2 (z)$ 2 + a(z)<& A + ■ ■ ■ need not concern 
us here, what is important is the second term describing the backreaction of the 
gauge fields on the brane geometry. In the antipodal Sakai-Sugimoto model C(z) = 
0 and there is no backreaction because of a geometric symmetry, but in the non- 
antipodal models C(z) ^ 0 and the scalar fields $, induced by the vector fields of 
the baryons, lead to attractive nuclear forced^H. The ratio of these attractive forces 
F a to the repulsive forces F r mediated by the abelian electric fields depends of the 
C(z) profile of the interaction term (fT71) . For the Sakai-Sugimoto models 


Fa 

F r 


= C 2 (z) 


1-C~ 3 /«o\ 
9 \u(z)J 


< 


< 1 , 


(18) 


b In our notations, the space-time indices 0,1,2,3, z of the effective five-dimensional theory are 
labeled = 0,1,2,3,21 while the space indices 1,2, 3, z of the same theory are labeled 

p, is ,... When we need the 9D indices for the whole D8 brane, we use m, n,... 
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so the net force is always repulsive. 

To see how that works, let us focus on the baryons in the Sakai-Sugimoto and 
other models of the holographic QCD. 


2.2. Baryons in holographic QCD models 

In the old hadronic string, the baryons were made out of Y-shaped configurations 
of three open strings connected to each other at one end; the other end of each 
string was connected to a quark. To realize this picture in holographic QCD, we 
need a baryon vertex (BV) — some object connected to N c open strings. The other 
ends of the strings must be connected to the flavor branes and act like the quarks; 
this would give the baryon its flavor quantum numbers. Witten had constructed 
the b aryo n vertex for the AdSs x S 5 model from a D5 brane wrapping the five- 
sphere-^ 41 generalized versions of this construction in Refs. [15] and [16] use Dp 
branes wrapping compact cycles carrying 0(N C ) Ramond-Ramond fluxes. 

In the Sakai-Sugimoto version of this construction, the baryon vertex is realized 
as a D4 brane wrapped around the S' 4 sphere (but localized in all other dimensions 
except the time). The S 4 carries N c units of the F 4 Ramond-Ramond flux, 


1 

(2tt) 3 Z| 



(19) 


so the Chern-Simons coupling of this flux to the U( 1) gauge field B living on the 
D4 brane acts a N c units of the net electric charge for the component B 3 : 


T 4 [ C 3 A e 2 ™' dB = N C [ B 0 dx° . (20) 

J D4 J 

In a compact space like S 4 , the net electric charge must vanish. To cancel the 
charge (®D we need to connect the D4 brane to open strings. The back end of an 
oriented open string has electric charge minus one, so we must connect the D4 brane 
with N c such strings; their front ends connect to the D8 flavor branes (since the 
strings do not have any other place to end) and act as N c quarks. 

We may put the D4 brane anywhere in space and anywhere on the cigar. How¬ 
ever, the S 4 volume, equivalently mass, increases with the u coordinate, so the 
lowest-energy location of the D4 is the cigar’s tip u = ua • At other locations, the 
brane feels a gravity-like force pulling it down to the tip. However, the strings con¬ 
nected to the BV pull it towards the flavor branes; in the non-antipodal models 
the D8 branes do not reach the cigar’s tip, so the strings pull the baryonic vertex 
up from the tip towards the lowest point uq of the flavor branes. The competition 
between the upward and downward forces on the BV determines its ultimate loca¬ 
tion. In some models, the forces reach equilibrium for the BV hanging on strings 
below the flavor branes^, while in many other models, including the non-antipodal 
Sakai-Sugimoto, the string forces win and the D4 sticks to the lowest point uq of 
the flavor branes^. 
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In all such models, the BV is a D p brane completely embedded in a stack of 
D(p + 4) flavor branes. As explained in Refs. [39l - l43l such D p branes are equivalent 
to zero-radius Yang-Mills instantons of the U(Nf) gauge symmetry on the flavor 
branes, and for Nf > 1 it may be smoothly inflated to a finite-radius instanton. In 
p+ 5 dimensions of the flavor branes, this instanton is a fat Dp brane wrapping some 
compact cycle, but once we dimensionally reduce to five dimensions, the instanton 
becomes a finite-size particle. Thus, in the low-energy effective five-dimensional 
theory of the holographic QCD, a baryon is realized as a finite-size instanton of the 
U(Nf) gauge theory. 

In the BV picture, each of the N c strings connecting the vertex to the flavor 
branes has electric charge 1 /Nf under the abelian U{ 1) subgroup of the U(Nf), 
so the whole baryon has abelian charge N c /Nf. In the instanton picture, the same 
electric charge is obtained from the CS coupling between the abelian electric field 
and the non-abelian magnetic fields of the instanton. 


Scs 


N c 

24tt 2 


tr 


( AT' 


2 - 4A 3 JF - -i-M 5 


10 


y)/ iA 

^ J d 5 x Aq(x) x I(x ), (21) 


where I{x) is the instanton number density of the magnetic fields (fTfil) . For Nf > 3, 
the CS couplings also endow instantons with non-abelian electric charges. 

The like-sign electric charges, both abelian and non-abelian, repel each other. 
It is the Coulomb repulsion between different parts of the same instanton that 
prevents it from collapsing to a point-like D-brane. However, the instantons do not 
grow large because the five-dimensional gauge coupling m decreases away from 
the z = 0 hyperplane: a large instanton would spread into regions of space where 
the coupling is weaker, and that would increase the instanton’s energy. Instead, the 
equilibrium radius of the instanton scales like 


t^inst 


1 

VXMkk 


( 22 ) 


For a holographic model of a baryon, this radius is unrealistically small. Indeed, 
using the p meson’s mass as a unit, the real-life baryon radius Rb ~ 3.4M“ , while 


in holography a <C M 1 ~ M KK . Moreover it raises the question of whether we may 


adequately describe such a small instanton using the DBI 
higher-order stringy corrections need to be included El 


CS action, or perhaps 


c In string theory, the DBI + CS action for the gauge fields on a D-brane is exact for constant 
tension fields Tmm however strong. But for the variable tension fields, the DBI action includes 
all powers of the Tmn but neglects their derivatives D^Tmn-, D^DpTmn-, etc. It is not clear what 
effect (if any) such higher-derivative terms would have on a small instanton. In a supersymmetric 
background, the instanton is BPS and its net mass is protected against stringy corrections, so the 
DBI action — or even the Yang-Mills action — gives the exact value. But what happens to small 
instanton in non-supersymmetric backgrounds is an open question. 
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On the other hand, assuming the DBI + CS description is valid, the small ra¬ 
dius (l 22 l) of the instanton allows for consistent expansion of the instanton related 
quantities in powers of 1/A.^° In particular, the leading contribution to the in- 
stanton’s mass M/ is 0(XN c Mkk) while the corrections due to ^-dependent gauge 
coupling and due to Coulomb self-repulsion are both 0(N c Mkk) ~ Mj/A. 

To see how that works, consider a static instanton — a time-independent con¬ 
figuration of SU(Nf) magnetic fields, plus the electric fields induced by the CS 
couplings, and the scalar fields induced by the tr (<&FF) coupling to the magnetic 
fields. Since the canonically normalized couplings in five dimensions are 0( A -1 / 2 ), 
the leading contribution to the instanton’s energy comes from a purely magnetic 
field configuration. In the DBI approximation 

Edbi = Jd 3 x.Jdz 2 ^ 2 1 ^ Str det [K{z)5 tlv + F^) - K 2 (z)^ , (23) 

where K(z ) = (27r a ') -1 gn{z) ~ 1 / {XM^ K £ A S ). For both 3 ym(^) and K(z) the de¬ 
pendence on the z coordinate is important only on the distance scale of order 1 /Mkkj 
so for instantons of much smaller size we may start with the approximation of con¬ 
stant K and constant five-dimensional coupling. In this approximation, the DBI 
energy is minimized by the magnetic fields that are exactly self-dual (with respect 
to the four-dimensional space of (a; 1 , x 2 , x 3 , z)); moreover, the DBI energy of an 
instanton is equal to its Yang-Mills energy: 

Edbi (instanton) = Eym (instanton) = — 5 —--— , (24) 

9ym\ z = 0 ) 

regardless of its radius and of the K parameter of the DBI action. 

The equilibrium radius of an instanton is determined by its (self-) interaction 
energy, in the next order in 1/A expansion. To this order, we assume the magnetic 
fields to be exactly self-dual — which allows us to use the YM action instead of DBI 
but the gauge coupling is ^-dependent, and we also account for the electric and 
the scalar fields. For small instantons, we may approximate the five-dimensional 
gauge coupling as 

= N c XM kk (B + DM^ k z 2 + 0(M* k z*)) (25) 

for some numerical constants B and D. For the Sakai-Sugimoto model B = £/277r 
and D = ( 8£ 3 — 5)/9£ 2 , while other models may have different values. Consequently, 
the YM energy of the non-abelian magnetic fields evaluates to 

-Ena = XN c Mkk ^B DM^^Z 2 + (26) 

where a is the instanton’s radius and Z is is the z coordinate of its center. 

The electric potentials Aq couple to the F llv F 111 ' products of the magnetic fields 
while the scalar potentials <h a couple to the F flu F llu . For the self-dual magnetic 


8tt 2 

Sym ( z ) 
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fields, both potentials couple to the same source J(x, z)Q“ v with Q “ l — the non- 
abelian electric charge, and the only difference is the coupling strength ratio C(z) 
( cf . Eqs. (flTl) and (fl8l) b For small instantons we may neglect the ^-dependence of 
this ratio and let C{z) ~ (7(0) = C, and as long as we do not go very far from 
the instanton (for distances r -C 1 /Mkk), we may also neglect the 2 -dependence 
of the gauge coupling. Consequently, both electric and scalar potentials become the 
four-dimensional Coulomb potentials 


$(x, z) = C Aq (x, 2 ) = 


2<3el 


BN c XMkk 


d 3 x' / dz' 


/(x', 2 ') 


(x-xO 2 + (*-*') 2 


( 27 > 


The electric fields lead to repulsive forces between the charges while the scalar 
forces lead to the attractive forces. Altogether, the net Coulomb energy amounts to 

B 


E C = 


N c XMkk j d 3 xdz tr((VA 0 ) 2 - (V$) 2 ) 

, l(x,z)xl(x',z') 


87 r 2 
(1 - C 2 )N t 


(28) 


d 3 x d 3 x' dz dz 1 


BXM J (x - x ') 2 + (z - z') 2 ’ 

and for a single instanton of radius a it evaluates to 

(1 - C 2 )N c 


E C = 


bBXMyj^a 2 


(29) 


Let us note that for C < 1 the electric fields are stronger than the scalar fields 
and the net Coulomb energy of the instanton is positive — which makes for a net 
self-repulsive force that prevents the instanton from shrinking to zero radius. In 
models with C > 1 (assuming they exist), the scalar fields would be stronger, the 
net Coulomb energy would be negative, which means a self-attractive force rather 
than self-repulsive. In such a model, the instanton would shrink to zero radius and 
our approximations would not be valid. 

Collecting the non-abelian (12^1) and the Coulomb (EHll self-interaction energies 
together we find the net energy of an instanton to be 

^T-inst = XN C M KK B + DM^ K Z 2 + DM^ K — + 5 _g^ 2 _^ 2 ^ a 2 ^ ■ 


Indeed, for a 2 ~ 1/(AM£ k ), both radius-dependent terms here are 0(1/A) correc¬ 
tions to the leading term. Minimizing the net energy, we find the equilibrium value 
of the instanton radius 


1 / 2(1 — C 2 ) \ 1/4 

AIkka/A \ 5 BD ) 


(31) 


The instanton center is in equilibrium at Z = 0 (the bottom of the U-shaped flavor 
branes), and the instanton’s mass is 


All—inst = N C M K K X B + 


2£>(1 - C 2 ) 


5 B 


0(1/A) 


(32) 




















16 V. Kaplunovsky, D. Melnikov, J. Sonnenschein 


For the antipodal Sakai-Sugimoto model 


a 


(1627T/5) 1 / 4 
Mkk \/A 


and 


M i_inst — 


N c M KK 



*JW + o( 1/A)) 


(33) 


while other models should have different 0(1) numeric factors. 

We may absorb two out of three model-dependent parameters B , C, D into a 
redefinition of the A and Mkk parameters of the effective five-dimensional theory, 
for example 

A —> ^ M = M K k \j^ , (34) 


thus 


8tt 2 


9ym( z ) 


= N c \m(i + M 2 z 2 + •••). 


(35) 


For static instanton or multi-instanton systems we may also get rid of the third 
model-dependent parameter C. Indeed, for the static systems $(x, z) = CA${x.,z), 
so the only effect of the scalar fields is to reduce the net Coulomb force by a constant 
factor (1 — C 2 ). We may simulate this effect without any scalar fields by using 
different gauge couplings for the electric and magnetic fields, g 2 { = (1 — C 2 )^ ag , 
or equivalently by rescaling the time dimension x° relative to the space dimensions 
(x 1 , x 2 , x 3 , z), 

t 


t 


but 


VI -c 2 

Consequently, the static instanton’s energy becomes 


(36) 


-Enet = N C \M + N C XM 3 J d 3 yidz /(x, z) x z 2 

+ -M I iz AT iJ * '(*'■*') + subleading (37) 

4AM J (x — x') 2 + {z- z ') 2 


= N C XM + N C \M 3 


rr2 a 

z + y 


N P . 


5A Ma 2 


+ subleading, 


hence in equilibrium Z = 0 and 


(2/5) 1 / 4 

mVa 


(38) 


Besides the radius a and Z , the instanton has other moduli the X 1 ’ 2 ’ 3 co¬ 
ordinates of the center (which corresponds to the baryon’s coordinates in three 
dimensions) and 4 Nf — 5 orientation moduli in the SU(Nf) gauge algebra. The net 
energy is degenerate in these moduli to all orders in 1/A. For finite N c — even if 
it’s very large — one should quantize the motion of the instanton in those moduli 
directions. Consequently, a holographic baryon acquires definite spin J and isospin 
/ quantum numbers: for Nf = 2 the baryons have / = J = 0,1,2,...,^ for even 
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Fig. 2. Conjectured phase diagram of large N c QCD.' — 

N c or 7=7=1,!,!,...,^ for odd i\r p5 | 24 |^ anc j there are similar spin-to-flavor 
multiplet relations for Nf > 2. 

However, in multi-baryon systems interactions between the baryons break the 
rotational and flavor symmetries of individual baryons, and only the overall SO( 3) 
and SU(Nf) symmetries remain unbroken. In holography, the multi-instanton sys¬ 
tems suffer from the same problem: the magnetic fields of multiple instantons in¬ 
terfere with each other, which spoils the degeneracy of the net energy with respect 
to orientations moduli of the individual instantons. In the large N c limit, this effect 
becomes more important than the quantum motion in the moduli space. 

Consequently, in this review we stick to classical static instantons with definite 
classical orientations in space and in SU(Nf). From the quantum point of view, such 
instantons are superpositions of states with different spins and isospins (or rather 
5'0(4) and SU(Nf) quantum numbers). At this level there is no difference between 
instantons being bosons or fermions. Since minimizing the classical energy of a 
classical multi-instanton system with respect to classical positions and orientations 
of all the instantons is already a very hard problem, going into the details of the 
quantization of moduli is beyond the scope of the review. 

3. Holographic Baryons at High Densities 
3.1. Nuclear matter in the large N c limit 

Let us first focus on the’t Hooft N c —> oc limit keeping the number of quark flavors 
small, Nf/N c <C 1. In such a case QCD perturbation theory is dominated by the 
planar gluon diagrams while contributions of the non-planar diagrams and of the 
quark loops are suppressed. In such a regime the T — g phase diagram is believed 
to look like Fig. [2] (For a review see Ref. (T2land references therein, also Ref. l44l ) 

Since the dynamics of the theory is dominated by the gluons, the quarks are 
sensitive to the gluonic background, but the backreaction from them to the gluons 
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is suppressed by Nf/N c . Consequently, at lower temperatures there is confinement, 
but for increasing temperature there is a first order transition to the deconfined 
phase. The transition temperature Ty is at the Aqcd scale and it is almost inde¬ 
pendent of the quark chemical potential fi q = /j,b/N c (as long as does not grow 
with N c ). 

It is not clear whether for N c —> oo the deconfining phase transition coincides 
with the chiral symmetry restoration for the light quarks. Several field theory ar¬ 
guments, for example in Refs 05} 06] and tHJ suggest that for /x q = 0 the two transi¬ 
tions should happen at the same point. However, these arguments do not work for 
fj-q > Opl and there are other arguments for the existence of confined but chirally 
restored phases, or deconfined phases where the chiral symmetry remains broken. 
In particular, some holographic models, e.g. in Ref. 1211 have deconfined but chirally 
broken phases even at zero chemical potential. 

For temperatures below the deconfining transition Ty and baryon chemical po¬ 
tentials below the baryon mass, equivalently for /i g < m q = Mt/N c , the thermal 
state of the theory is a dilute gas of glueballs and mesons, whose interactions are 
suppressed by powers of N Cl and almost no baryons or antibaryons: the masses of 
mesons and glueballs are O(Aqcd), and they scale as Af?; on the other hand, the 
baryons, which are made out of N c quarks, have masses of order -/V c Aqcd, so their 
relative abundance in thermal equilibrium is exponentially suppressed. The interac¬ 
tion energy between baryons also scales as N c . Hence at fi q ss m q there is an abrupt 
phase transition to the bulk nuclear matter with finite baryon density. 

Unlike the ordinary nuclear matter, which is in the quantum liquid state, the 
large N c nuclear matter is crystalline solid since the ratio of kinetic energy to poten¬ 
tial energy decreases with N c . Indeed, the potential energy of baryon-baryon forces 
scales like N c ; more precisely,^ in the large N c limit the two-baryon potential 
becomes 


V ~ N c x Ac(r) + N c x A s (r) (U • I 2 ) (Ji ■ J 2 ) + 

+ N c x A T (r) (Ir • I 2 ) (3 (n • Jr) (n • J 2 ) - (Ji • J 2 )) + O (1 /N c ) , (39) 

for some ^-independent profiles Ac(r), As(r), and At(t) for the central, spin- 
spin, and tensor forces; their overall magnitudes are A ~ Aqcd for r ~ 1/Aqcd- 
Classically, this potential tries to organize the baryons into some kind of a crystal, 
where the distances between neighboring baryons do not depend on the N c , while 
the binding energy (per baryon) scales like -/V c Aqcd- Quantum-mechanically, the 
baryons in such a crystal behave like atoms in ordinary crystals: they oscillate in 
their potential wells with zero-point kinetic energies 


K 


7T 

2m b d 2 


Aqcd 
N c d 2 ’ 


(40) 


where d ~ 1/Aqcd is the ^-independent diameter of the potential well. Therefore, 
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at zero temperature the ratio of kinetic energy to the potential energy scales like 


I< 1 
V ~ N* 


(41) 


and becomes very small for large N c . 

At higher temperatures the kinetic energies of baryons become larger, K ~ T, 
but in the confined phase we are limited to T < Td ~ Aqcd- Consequently, the 
kinetic to potential energy ratio scales at most like 


K 1 

v ~ T c ’ 


(42) 


which is larger than (TUI) but still small in the large N c limit. Consequently, for 
large N c neither zero-point quantum motion, nor thermal motion of baryons can 
destroy the baryon crystal, so the nuclear matter remains solid all the way to the 
deconhning temperature. 

Before holography, the best models for the large-iV c nuclear crystals were lat¬ 
tices of skyrmionsl- In this framework, Goldhaber and Manton had found a curious 
phase transition from a lattice of whole skyrmion s at low chemical potential to a 
denser lattice of half-skyrmions at higher potential.^ (According to Refs l9l and [TUI 
the whole-skyrmion lattice at low potentials has face-centered cubic (fee) arrange¬ 
ment while the half-skyrmion lattice at higher potential is simple cubic.) In the 
half-skyrmion-lattice phase, the order parameter for the chiral symmetry breaking 
vanishes after space averaging, so in QCD terms the transition is interpreted as 
chiral symmetry restoration at high fi q . 

In QCD with N c = 3 and two massless flavors there is a similar chirally- 
symmetric phase at low T and high fi q . This phase is a quark liquid rather than a 
baryon liquid (the quarks are no longer confined to individual baryons) and there 
is a condensate of quark pairs making this liquid a color superconductor. But for 
large N c the situation is more complicated: there is no color superconductivity, and 
there is no deconfinement for T < Ty. Instead, the dense cold nuclear matter forms 
a phase which combines the features of the baryonic and quark phases: the quarks 
fill up a Fermi sea, but the interactions near the Fermi surface are strong, so the 
excitations are not free quarks or holes but rather meson-like quark hole pairs or 
baryon-like states of N c quarks. 

For n q Aqcd, the interior of the Fermi see is chirally symmetric, but near 
the Fermi surface the symmetry is broken by the chiral density waves.® (Although 
Son and Shustei®] argue that such waves develop only for very large N c > 10 3 .) 
To be precise, the chiral density waves mix the chiral symmetry of the quark phase 
with the translational symmetry rather than simply break it. Averaging over space 
restores the chiral symmetry, just like it happens for the lattice of half-skyrmions. 

On the other hand, for jj q just above m q = My/N c , the baryonic crystal has 
a completely broken chiral symmetry. Thus, at some critical /iq C> = 0(m q ), there 
should be a chiral symmetry restoring phase transition from the baryonic crystal 
to a distinct quark phase. In holography, this transition is believed to be dual to 
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Fig. 3. Two-body nuclear potential in holographic QCD. 

the “popcorn transition” from a three dimensional to a four dimensional instanton 
lattice.™ 


3.2. Effect of the large A limit on the holographic nuclear matter 

In holography, the semiclassical description of the gravity side of the gauge/gravity 
duality in terms of metric, fluxes, and branes requires the limits of large N c and 
also large’t Hooft coupling A = gy M N c . In the large A limit, the baryons become 
very heavy: in units of the mesonic mass scale (M ~ Aqcd) the baryon mass is 
Mb ~ A N C M. However, the interactions between the baryons do not grow with A: 
even for two baryons right on top of each other, the repulsive potential between them 
is only V ~ N C M ~ Mb/ A. At larger distances, the forces are even weaker since the 
hard-core radius of a holographic baryon shrinks with A as Rb ~ M -1 A -1 / 2 . Outside 
this radius, the repulsive potential decreases as 1/r 2 until r ~ M , at which point 
it becomes dependent on the meson mass spectrum of a specific holographic model: 
In some models, the potential becomes attractive for r > M -1 while in others it 
remains repulsive at all distances. The overall picture is shown in figure H3 

Since the nuclear forces are so weak in the holographic QCD, all transitions 
between different phases of the cold nuclear matter happen at chemical potentials 
p very close to the baryon mass: just below Mb we have glueball/meson gas (or 
vacuum for T = 0), while just above Mb we have dense quark matter. To see the 
baryonic matter phase (or any other intermediate phases) we need to zoom into the 
p ~ Mb region. 

Figure |4] illustrates this point: to see the baryonic matter phase between the 
vacuum and the quark-matter phases on the plot of baryon density p as a function of 
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Fig. 4. Density as a function of the chemical potential in the large A regime for attracting (left) 
and repelling (right) baryons. The nuclear matter phase is confined to a narrow window of the 
order A// ~ MfJ A. In the naive diagram the transition occurs directly from the no-baryon to a 
quark phase. 

the chemical potential /x, we need to zoom into the narrow range ft—My = O^AIy/X). 
The figure also shows that the thermodynamic order of the phase transition between 
the vacuum (or the meson/glueball gas for T > 0) and the baryonic matter depends 
on the sign of the long-distances nuclear force. If the force becomes attractive at 
long distances, then bulk baryonic matter exists at zero external pressure and has 
fi > Mb — -Ebind, where -Ehind is the binding energy. Consequently, the transition 
from the vacuum (or gas) to the nuclear matter (in the form of a baryonic crystal) 
is first order as shown on figure 0J left). On the other hand, if the nuclear forces are 
repulsive at all distances, then the bulk nuclear matter does not exists except at 
positive external pressures, and its chemical potential must be /x > Mb. Moreover, 
/x raises monotonically with the pressure and the density, so the transition from the 
vacuum to the bulk nuclear matter is second order as shown on figure IDd'ight) ■ 
Phase diagram of the Sakai-Sugimoto model was studied, for example, in 
Refs. [50H521 and jj?2J where the narrow nuclear matter phase was ignored. In the 
latter examples the scenario on the right part of figure 0 ] applies. Here we make no 
assumptions about the long-distance nuclear forces. Consequently, we cannot say 
anything specific about the transition from the vacuum to the nuclear matter phase. 
Instead, we focus on the transition from that phase to the quark liquid phase, which 
correspond in the holographic picture to changing the configuration of the instan¬ 
ton crystal, from a three-dimensional lattice, through a sequence of intermediate 
steps, to a four-dimensional one. Or rather, in this article, we review a simpli¬ 
fied problem, namely the transitions between one-dimensional and two-dimensional 
instanton lattices. 


3.3. Multi B ary on Systems in holographic QCD 

In the large N c limit, nuclear forces between the baryons are dominated by the 
static potentials. Holographically, a static system of A baryons corresponds to a 
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time-independent configuration of the non-abelian magnetic flavor fields F° v (x.,z) 
(/x, v = 1, 2, 3, z) of net instanton number A, 



167T 2 


tr {FkxF^) 


= A, 


(43) 


accompanied by the Coulomb electric Ag(x, z) and scalar <3> a (x,z) potentials in¬ 
duced by their Chern-Simons and couplings to the magnetic fields (see sec¬ 
tion [2T2] and Ref. [28] for the details). The whole configuration should minimize the 
net DBI + CS energy of the system subject to the constraint (l43l) . 

In the A —> 00 limit, the DBI energy of the magnetic fields is 0{ A) while the net 
effect of the electric and scalar fields is only 0(1). Moreover, the magnetic fields are 
concentrated within 0(A -1 / 2 ) distance from the 2 = 0 hyperplane, so to the leading 
order we may approximate the 4 + 1 dimensional spacetime as flat. Thus Eq. (1231) 
becomes 




Edbi 


9ym 


d 3 xdz Str 


\J det(A'o<5 flU + Ff, „) — I \0 


(44) 


where we neglect the 2 dependence of the gy M and Ko = K(0). Similar to the Yang- 
Mills energy of an A-instanton system, this leading-order DBI energy is minimized 
by the self-dual configurations of the magnetic fields (x, 2 ). In fact, all such self¬ 
dual configurations (of the same instanton number A) have the same leading-order 
energies 


-Elo = = 5 — x A = AXN C M , (45) 

3ym 

and Ad was introduced in (1351) . The self-dual configurations form a continuous fam¬ 
ily parameterized by 4iV/A moduli, which correspond to the locations, radii, and 
SU(Nf) orientations of the A instantons. But the leading-order energy (l45l) does 
not depend on any of these moduli. 

Fortunately, the sub-leading corrections to the net energy lift the degeneracy of 
the leading order, which provides for the 0(A°) interactions between the baryons. 
To work out such interactions we need the degenerate perturbation theory for the 
magnetic field configurations and their energies. Let us outline the formal procedure. 

At first order of the perturbation theory, one (a) limits the F“^(x, 2 ) configu¬ 
rations to the degenerate minima of the leading-order energy function, i.e. to the 
self-dual magnetic fields; (b) calculates the 0(A°) corrections A E for the energies 
of these configurations; (c) minimizes A E among the self-dual configurations. At 
the next order, to compute the ©(A” 1 ) corrections to the energy, one needs to find 
the 0(A -1 ) corrections AF“„(x, 2 ) to the magnetic fields, which would no longer 
be self-dual. This is however beyond the aim of this review. 

In more detail, to get 0(A°) interactions between A baryons, we proceed as 
follows: 


(1) Use Atiyah-Drinfeld -Hitchin-Manin (ADHM) construction^ to obtain gen¬ 
eral self-dual magnetic field configurations, for a review see e.g. Ref. [54] Those 
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are encoded in terms of A x A and A x Nf matrices obeying certain quadratic 

(ADHM) constraints. The first task is to solve these constraints and write down 

the ADHM matrices in terms of the instantons’ locations, radii, and orienta¬ 
tions. 

Given the ADHM matrices, work out the instanton number density profile 

I (x,z) = -^2 tr(F KA F M1/ ), (46) 

for Nf > 2 one would also need the non-abelian adjoint density 

Hx, x) = ^d abc F b x F^. (47) 

Next, calculate the O(X 0 N c ) corrections to the net energy of the system, which 
include the effect of the ^-dependence of the five-dimensional gauge coupling 
and the Coulomb electric and scalar potentials induced by the Chern-Simons 
and the &FF couplings, thus 

A E = A Ena + A Eq ■ (48) 


The ^-dependent gauge coupling changes the DBI energy of the magnetic fields 
by 

AAna = N C M f d 3 xdz XM 2 z 2 x J(x, 2 ), (49) 


while the Coulomb energy depends on the Nf. For Nf = 2, the U( 2) CS and 
&FF terms couple the SU(2) magnetic fields to the 17(1) Coulomb fields only. 
Consequently, the Ag and the <f> a fields are abelian and couple to the instanton 
density J(x, z). Their net energy is simply 4+1 dimensional Coulomb energy 


A Ec 



d 3 x 2 dz 2 


7(xi,zi) x J(x 2 ,z 2 ) 
AM 2 ((xi - x 2 ) 2 + (zi - z 2 ) 2 ) ' 


(50) 


For Nf > 2, the U(Np) CS and terms couple the SU(Np) magnetic fields 
to both abelian and non-abelian electric and scalar fields; the abelian Coulomb 
fields are sourced by the instanton density I (x, z) while the non-abelian fields 
are sourced by the adjoint density 7 a (x, z). Altogether, the net energy of these 
Coulomb fields is 


A Ec 



2/(xi, zi) x /(x 2 ,2 2 ) 
AA/M 2 ((xi - x 2 ) 2 + (21 - z 2 ) 2 ) 
4J a (xi,zi) x J a (x 2 ,z 2 ) ) 

+ AM 2 (( Xl - x 2 ) 2 + ( 2l - z 2 ) 2 ) j ' j 



The previous three steps give A E as a function of baryons’ locations, radii, 
and SU ( Nf ) orientations. In the final step minimize the A E with respect to all 
these moduli. 



24 V. Kaplunovsky, D. Melnikov, J. Sonnenschein 


This four-step procedure is fairly straightforward for a few baryons — cf. cal- 
dilations of the two-body nuclear forces by Kim and Zahed^l and by Hash imoto 
et a/t'®]However, it becomes prohibitively difficult for large numbers of baryoni^ 53 
and outright impossible for infinite baryon crystals. At best, one can survey a small 
subspace of the A-instanton moduli space and try to minimize the A E over that 
subspace. For example, we may assume that all the instantons have the same ra¬ 
dius, that their centers form a periodic lattice of some particular symmetry, and 
that the orientations of the instantons also form some kind of a periodic pattern. 
This gives us an ansatz for all the A x 4Nf moduli in terms of just a few overall 
parameters, and we can try to calculate and minimize the A E as a function of 
these parameters. However, any such ansatz is likely to miss the true lowest-energy 
configuration of the system. Indeed, in condensed matter guessing the crystalline 
symmetry of some substance from the properties of the individual atoms is a game 
of chance with poor odds, and there is no reason why the instanton crystals should 
be any simpler. Moreover, even if we could somehow guess all the symmetries of the 
instanton crystal, actually working through above four steps is impossible without 
additional approximations (besides A > 1). 

In the examples we consider in the remainder of this review we will use several 
approximations. In section 13.41 we will illustrate our main expectations about the 
high density properties of the holographic nuclear matter by considering point-like 
instantons. In section ITTl where the low-density will be presented, we will need the 
sparse lattice approximation (the size of the instanton is much smaller than the unit 
cell size) in order to compute the Coulomb energy. In section I4~21 this approximation 
will be necessary to solve for the “high-density” phase. In fact, this is as far as the 
analytical methods have allowed us so far.In section[5]we will further neglect the 
contribution of the many-body forces to the instantons’ interactions. This allows for 
an efficient numerical analysis and has lead to a plethora of high-density phases!^ 
In both cases, however, we will restrict only to the one-dimensional crystals. 

To make the baryons form a one-dimensional lattice instead of spreading out in 
three dimensions, we curve the two transverse dimensions of the flavor branes, say 
X 2 and X 3 , similar to the curvature of the 2 coordinate. In terms of the effective 4+1 
dimensional theory, this corresponds to the five-dimensional flavor gauge coupling 
depending on the X 2 and the X 3 as well as the X 4 = z, 

o_2 , v 

2 1 = XN C M[1 + M%x\ + M\x\ + Mlxl + 0(l¥) . (52) 

This gauge coupling acts as a harmonic potential for the instantons which pulls 
them towards the X\ axis, so at low densities the instantons form a one-dimensional 
lattice along the x\. At higher densities, the instantons push each other away from 
the X\ axis and form more complicated two-dimensional or three-dimensional lat¬ 
tices, starting with the zigzag, the configuration shown in figure [5] To make sure 
the transition from the straight chain to the zigzag happens for lattice spacings 
much larger than the instanton radius (which is required by the two-body force 
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Fig. 5. The zigzag configuration is the leading order small deformation of the one-dimensional 
chain. 


approximation), we will assume M 2 <C M 3 , M 4 . 


3.4. Point charge approximation 

In order to illustrate our expectation about the phase structure of holographic 
baryons at high density let us consider a drastic simplification. Let us take the limit 
in which the size of instantons goes to zero, in which they become like point charges. 
In this limit the only surviving moduli of the instantons are their positions. As 
explained in the previous section, by tuning a potential in the transverse directions 
X 2 and £ 3 , see Eq. (ED, we will confine the lattice of instantons-charges to a one¬ 
dimensional chain along x\. We will assume that potential in one of the directions, 
say X 2 , is weaker than in the others, specifically M 2 <C M 3 l M±. Later, in section[D 
we will observe that keeping the ration M 3 /M 4 non-trivial allows for a number of 
new high-density phases. 

In the point-charge limit there are two competing forces exerted on instantons. 
One comes from the curvature-generated “non-abelian” potential, which confines 
them to the x\ dimension, and the second is the four-dimensional Coulomb re¬ 
pulsion, which dislikes instantons sitting too close to each other. At low densities 
the equilibrium configuration must be a straight periodic chain of instantons. The 
one-dimensional density is 1/D, where D is the lattice spacing. 

When D goes to zero, equivalently density grows high, the Coulomb repulsion 
will force the instantons to leave the one-dimensional alignment and to expand in 
the transverse directions at the expense of increasing their non abelian energy. The 
leading displacement mode is the zigzag shown in figure [U in which each pair of 
nearest neighbors goes in the opposite directions away from the chain. Recall that 
from the point of view of the flavor branes the neighboring instantons in the zigzag 
move to the opposite branches of the U-shaped configuration. This was called the 
popcorn transition in Ref. 1291 

Let us now study this transition quantitatively. We replace the instanton density 
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I ( x ) by the sum of delta-functions 

OO 

'(*)= £ ^ 4 ) (x-nD), (53) 

n =—oo 

where D is a 4-vector generating translations from one site of the chain to a neigh¬ 
boring one, here chosen to be along x\. For the straight chain the non-abelian part 
of the energy, calculated per instanton, gives 


d 

-Ena = N C \M J dx i J d 3 x I(x) (l + M 2 a ; 2 + M|a ; 2 + M 2 ;c 2 ) = 

o 

= N C XM (1 + M%xl + Mlx\ + M 2 x \). (54) 

This energy is minimized for the choice X 2 = £3 = X 4 = 0. From Eq. m th e 
Coulomb energy per instanton is given by the sum 


_ N c yv 1 
“ 4AM ^ (nD ) 2 

njt 0 V ' 


N c 7T 2 
AM 12E 2 ' 


(55) 


Let us investigate stability of the charges in the straight chain against the tran¬ 
sition to the zigzag. The preferred zigzag deformation is in the X 2 direction as we 
have frozen X 3 = £4 = 0 by requiring M 2 ■'C For every charge in the 

zigzag A ^2 = ±e. As a result the averaged change in the non-abelian energy (per 
instanton) will be 


AE na = iV c AMM|e 2 . 


(56) 


It is also straightforward to evaluate the Coulomb energy per instanton in the zigzag 
phase: 


Ec = 


N c 

4AM 


E ■ 

1 even n^O 


iDy 


E 

odd n 


(nd ) 2 + (2e) 2 


N r . 


(- 


+ 


7T 7T6 A 

- tanh — 

16eE D J 


(57) 


AM \48-D 2 

For small e let us expand the energy difference between the zigzag and the 
straight phases: 

N r . 


A E = 7V c AMM|e 2 


7r 4 e 2 


7r 6 e 4 


AM V 48E 4 120 E 6 


0(e b ) 


(58) 


The sign of the e 2 -term depends on the density (lattice spacing). It is positive 
for small densities and negative for large. Therefore there is a second order phase 
transition, in which e acquires a non-trivial vev, i.e. the straight chain turns into a 
zigzag. The critical density corresponds to the point, where e 2 -term changes sign: 

7T 


D = D c = 


2 • 3 1 / 4 v / MM 2 A ' 


(59) 
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For the spacing slightly smaller than D c the order (zigzag) parameter has the mean 
field behavior: 



(60) 


More generally e is given by the solution to 


(ire/D) 3 cosh 2 (ne/D) 3 4 



(61) 


sinh(27re/-D) — (2ire/D) 4 \ D J 


The zigzag amplitude grows monotonously with the density. However at larger den¬ 
sities, that is smaller spacing between instantons in xi, it may become favorable for 
instantons to organize themselves in a more complicated configuration with respect 
to X 2 in order to minimize the Coulomb energy. Indeed, as was shown in Ref. 1291 
the zigzag first turns into a double zigzag (four layers in £ 3 ), then to a three-layer 
configuration etc. Gradually the original one-dimensional lattice turns into a two- 
dimensional one. Except for the first transition to the zigzag, which is second order, 
the remaining ones are of first order. 

A similar analysis can be done for two- and three-dimensional lattices. For the 
three-dimensional case one does not need additional parameters M2, M3, there is 
only the curvature of the holographic z-dimension, M 4 = M. If one considers a 
simple cubic (sc) lattice, it is likely that at large densities it will split along the 2 
direction into a three-dimensional analog of the zigzag. In such a transition every 
even site will displace in the ±7 direction while the every odd one will displace in 
the opposite. In this second order phase transition the original sc lattice will split 
into two fee layers. 

The transition sc —> 2fcc can be studied quantitatively. The only difference is 
that in two and three dimensions the Coulomb energy Eq per instanton will diverge. 
However, what matter for the the question of stability is the energy difference of 
the two configurations, which is finite. Overall, one finds 



(62) 


where the following dimensionless quantities have been introduced 



(63) 


This gives the critical density corresponding to the lattice spacing 



(64) 


The sc lattice however is unlikely the ground state of the system of point charges. 
The minimum energy configuration of in this case is unknown, but we believe it to 
close packing, that is the largest interatomic distance between the nearest neighbors 


for a given density, which is achieved in the fee lattice. Stability analysis of the fee 


lattice is trickier because for the fee lattice there is no natural way to split into two 
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sublattices. One can expect two things to happen when the density is increased. 
First there can be a first order transition to a multi-layer four-dimensional lattice. 
Second the splitting into two sublattices can occur through a breaking of the cubic 
symmetries (half of nearest neighbors will go one way and another half the other 
way). The latter will be restored as soon as the separation between the sublattices 
will be large enough. 

The example of the point charges demonstrates the main expectations about the 
phase structure of the holographic baryon crystal. When squeezed hard enough the 
lattices of baryons-instantons will expand in transverse dimensions to form higher¬ 
dimensional structures. In particular, the three-dimensional lattice of instantons 
will turn into a four-dimensional one with a non-zero thickness in the holographic 
dimension. As the density grows high the number of possible phases (lattice config¬ 
urations) grows rapidly. This number is even higher if one goes beyond the point- 
charge approximation, as one has to take into account other instantons’ moduli, 
such as orientation. In the remaining sections we will review the current results on 
the phase diagram of the nuclear matter phase in the generalized Sakai-Sugimoto 
model, so far obtained for the one-dimensional instanton lattices. 

4. Exact solutions 

In this section we will discuss some analytical Nf = 2 instanton solutions, which de¬ 
scribe holographic one-dimensional baryon crystals close to the threshold chemical 
potential, i. e. around the nuclear-to-quark matter phase transition. As we discussed 
in the previous section for analytical treatment to be viable, we need some assump¬ 
tions about the instanton moduli. In this section we will assume a certain periodicity 
of the instantons’ orientations. In order to show that the solutions, considered here, 
are true ground states it is desirable to expand the available moduli space. This will 
be done in section 0 where a larger class of solutions will be scanned numerically 
based on the further two-body force approximation. 

4.1. Straight chain 

Let us start from reviewing the straight chain of instantons, which is the appropriate 
solution for the low-density baryon crystal. The solution for an infinite periodic 
chain of SU( 2) instantons was first obtained by Harrington and Shepard in Ref. l59l 
That solution corresponded to parallel-oriented instantons. Two decades later, in 
Ref. l60l Kraan and van Baal derived a generalization of the Harrington-Shepard 
solution, which described an instanton chain with a constant relative orientation 
shift between any pair of neighboring instantons. In this subsection we will review 
the result of Ref. [60]in the context of holographic QCD. 

One reason why it took a while to generalize the solution of Ref. [59] to the 
solution of Ref. I60lwas that the former can be obtained through the simple ’t Hooft 
ansatz, while the case of instantons of variable orientation required the details of 
the ADHM construction. ADHM data are matrices, which contain the information 
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about the locations of the instanton centers, their radii and orientations. For an 
infinite chain the matrices are infinite-dimensional. Specifically, Sp(k) A-instanton 
solution is encoded by two quaternionic matrices: an A x A symmetric matrix X 
and a k x A vector Y, which satisfy the constraint 

x^x + y f y is real symmetric. (65) 


In the case of Sp(l) ~ SU{2) the solution can be reformulated in terms of four 
real symmetric matrices and real vectors Y M : 

x = = r 4 + ir j r j , y = y^ = y 4 + ir j Y j , (66) 


where r J , j = 1,2,3 are Pauli matrices. The real matrices satisfy the constraint 

f-flVK X 

[r^r"] + y »® y v - y v ^ y > 1 = —([r K ,r A ] + y k ® y x - y*«y K ), (67) 

where ( Y M ® Y v ) mn = YJ£ x Y". Physically, the diagonal matrix elements r^ n are 
the four-dimensional coordinates of the instantons’ centers. The components Y£ 
combine the radii and the SU( 2) orientations of the instantons. In particular, we 
introduce a n = \Y n \ for the real radii, and A SU{2) matrices y n (equivalent to 
unimodular quaternions Y n /a n ) parameterizing the instanton’s orientations. Equa¬ 
tions © can be equivalently written as 

^[T^rlrrm + S "*“" X tr (2/mJ/n(-”'“)) = 0, (68) 

where 77 “^ is the’t Hooft’s symbol mapping between the SU(2) gauge and the SU(2)l 
inside Spin(4) = SU(2)l x SU(2 )r, 


<1 = 0, < = -<, Vu = + 8 ?, Vij=£ aij , a,i,j = 1,2,3. 


(69) 


For a given ansatz, with selected radii, positions and orientations of instantons, 
the off-diagonal matrix elements r^ n = a^ n are determined through the ADHM 
conditions ©. The ADHM data are somewhat redundant — an O(A) symmetry 
acting on all the F^ n and Y n = a n y n does not change any physical properties 
of the multi-instanton data. This symmetry includes Z^ which flips the instanton 
orientations y n —y n (independently for each n). It also includes small 50(A) 
rotations that change the off-diagonal elements by Sal^ n = e mn (X m — X n ) M + 0(e 2 ). 
To eliminate these rotations, the ADHM equations (l68l) for the off-diagonal elements 
should be combined with additional constrains (one for each m ^ n), for example 


Vm^n: (X m - X n Ya^ n = 0. 


(70) 


To find the self-dual vector- potential one also needs to solve an additional set 
of matrix equations. Luckily, for the kind of analysis we do here what we need is 
just the expression for the instanton density. In terms of the matrices and Y M 
it can be constructed as follows. Provided that (l67l) is satisfied define a real Ax A 
symmetric matrix 

Lmn(x) = 'y ] (r% m — + — a m a n tr^y^y ^. (71) 

t 
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The instanton density is then given bjP 


I{x) = -Y^nniogdet(L(a;)). 


(72) 


We will look for a solution, which corresponds to an array of equal size instan- 
tons arranged in a straight periodic lattice. Periodicity means discrete translational 
symmetry S : aT 1 —> x 11 + D M of the ADHM solution; in the language of the 
and Y£, this symmetry acts as 


T* S = P M + 1 to keep the 1 — r^ 1 invariant 

(I'M ->• Y, G{ y^) S mn = (yy), 

m 


(73) 

(74) 


for some 0(A = oo) matrix S mn and SU(2) matrix G. 

Physically, G rotates the orientation of an instanton relative to its immediate 
neighbor. The more distant neighbors are related by G n ~ m rotations, which generate 
a 17(1) subgroup of the SU(2 ). Without loss of generality, we take G = exp (icjyrz/fy 
for some “twist” angle (f> between 0 and 27 t. 

In this section wc will take the direction of the instanton chain to be the x 4 , 
while the transverse directions are x 1 j x 2 , and x 3 . We will also make X 3 the preferred 
direction for the zigzag deformationo In terms of equation (El this means D M = 
(0,0,0,77). Finally, we take S mt „ = <5 m ,n+ 1 (shifts from the n — th instanton to 
the (n + 1) — st) and Yq = (0,0,0, a) (where a is the radius of the instantons). 
Consequently, the translational symmetry El El requires 

ir^Y£ = a exp <*=> Y£ = (0,0,asin(m/>/2),acos(n^>/2)) (75) 

and 


T mn = Dd^xnSmn + f M (m-n), (76) 

where the T M (m — n) do not have separate dependence on m and n but only on 
m — n. Combining these symmetry conditions with the ADHM constraint (1671) . we 
get 


Y£ ® Y. if* 


r 

r 


4 

mn 

1 

mn 


r 


3 

mn 


a 2 cos [(to — n)<t>/2 ], 
Dn 5 mn , 

r 2 =0 

mn 

a 2 sin [(to — n)<j>/2 ] 
d to — n 


for to ^ n, but 0 for to = n. 


(77) 

(78) 

(79) 

(80) 


d Notice that this is different from the convention used in the previous section, where the chain 
was along x 1 and the zigzag was formed in x 2 . We will return again to that convention in the next 
section. 
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To calculate the instanton density of the periodic chain we need the determinant 
of the infinite matrix L CD- This determinant is badly divergent, but we may obtain 
it up to an overall infinite-but-constant factor from the derivatives 


d/j, logdet(L(a;)) = 2tr ((/l — T M )L 1 (x)) . (81) 


For three of these derivatives (/i = 1,2,3) the trace converges, while for [i = 4 
the trace diverges, but can be regularized using symmetry X4 —> — X4, n —> —n. To 
evaluate those traces, it is natural to use the Fourier transform from infinite matrices 
to linear operators acting on periodic functions of 9 (mod 2tt)EI Consequently, L 
becomes 


L 

where T(Q) 
and r 3 ((9) 


x i + x \ + ( x 3 ~ T 3 (0)) + ^X4 + iD—' S j + T{9), 

na 2 S(9—^) + 7ra 2 (5(d+^), 

7ra 2 f 1 — ^ for 0 < 9 < ^ and 2n — ^ < 9 < 27r, 
~°l forf<fl<27r-f, 


(82) 

(83) 

(84) 


and L 1 becomes the Green’s function of this operator. Calculating this Green’s 
function one obtains the following expression upon integration of the traces (IHTl) : 


det(L) = (cosh^-!- + 


D 


7 rar . <f>ri 
—— sinh ——- 
Dr 1 D 


cosh 


where 


rf + r 2 ~ (7ra 2 /D) 2 

2nr 2 

2-KX4 

1 T ’ 


+ 


— cos ■ 


<t>ri 


sinh- sinh 

D 


(2n - (j))r 2 _ 
D 

(27r - (j))r 2 


D 


na" . , (27t — (f>)r 2 


Dr 2 


sinh ■ 


D 


(85) 


r = Xi 


r 0 = x. 




( , a2 (^ 

+ x 2 

+ 

[ X:i + 2 



( a 2 cj)\ 

+ x 2 

+ 

l3+ 2 D, 


( 86 ) 


This is precisely the result obtain by Kraan and van Baal in Ref. [501 . 

Let us now put this instanton solution in the context of holographic baryons. 
In line with the program discussed in section [3] we are going to evaluate 0(A°) 
corrections to the leading order energy result. For that we will plug the density of 
the “free” instanton solution into equations (l49l) and (l50l) . Then we will minimize 
the energy with respect to the relevant moduli, which in this case are the radius 
and the twist angle. 

The expression for the instanton density follows from (l72l) after plugging so¬ 
lution (l85l) . Although the expression derived from (l85l) is way too complicated to 


'This operation is also known as Nahm transform, cf. Ref. l62l 
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print, several moments of the instanton density may be obtained via integrating by 
parts, thanks to the double D’Alembertian in Eq. d): 


d 4 " 


~xl(x) = A, 

J d 4 xl{x) x x v = {x v ) = tr(r y ), 

J d 4 xl(x) Xi¥ =(x^x v )= tr(r^r^) + id^tr(T) 


In particular, 


where T mn = - a m a n tr(y^y n ). 


<*?> = <* 2 > = 


— exactly as for a single stand-alone instanton of radius a, but 


<*§> = y + 


AD 2 


(j)( 2 n - 0), 


(87) 

( 88 ) 

(89) 

(90) 

(91) 

(92) 


where the second term is due to interference between the instantons. Curiously, 
the interference term vanishes for <j> = 0, i .e. for instantons with the same SU(2) 
orientations. 

The calculated moments of the instanton density yield the leading corrections 
to the non-abelian part of the energy (PI) . Taking the expansion of the five¬ 
dimensional gauge coupling from Eq. (l52l) and using Eqs. (ED and (l92l) . assuming 
small size instantons (a -C M _1 ), one obtains, 


AEna = N C XM 


(M- 


Ml 


Ml) 


Mia 4 . 


(93) 


Here M\ = M 4 is appears since we stretch the chain along the X 4 dimension. Later 
in this section we will assume Mi = M 2 = M' and M 3 = M. 

Adapting Eq. (15(71) one finds the abelian Coulomb energy per instanton to be 

Ec = 256n 2 X M l dx4jd 3 x(d,Dlogdet(L)) 2 . (94) 

For generic lattice spacings D , this integral is too hard to take analytically. But it 
becomes much simpler in the D a limit of well-separated instantons. 

For large lattice spacing D a the leading correction to the Coulomb en¬ 
ergy EH) evaluates to 

1 | 4-7T 2 + 3(7T — (j)) 2 

5 a 2 + 


AE ° “ H 


30.D 2 


+ 0 {a 2 /D 4 ) 


(95) 


Now we need to minimize the total energy with respect to the moduli: the 
instantons’ radius a and the twist angle (j). Combining (l93l) and (l95l) . we find the 
minimum at 


= a 0 - 


n 2 <4 

12 D 2 


+ O(a 4 0 /D 2 ), 


'min 


(96) 
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where ao is the equilibrium radius of a standalone instanton (1381) . We see that 
in the “low-density” regime the instantons in the straight chain prefer the anti¬ 
ferromagnetic order, i.e. any two neighboring instantons have an opposite orienta¬ 
tion. 


4.2. Zigzag 

Let us now address the question of the stability (or instability) of the straight 
instanton chain considered in the previous subsection against transverse motion (in 
34 , 2,3 directions) of the instantons that would get them out of the linear alignment. 
If such a motion can decrease the chain’s net energy, than the chain is unstable. 

Following the conventions of this section, we set M\ = M 2 = Ad' and M 3 = M 
and assume M < Ad', so that the leading instability would be in the £3 direc¬ 
tion. The anisotropy also breaks the degeneracy between different directions n of 
the 517(2) twist exp(i<^(n • r)/2) between adjacent instantons. The lowest-energy 
direction of the twist is now T 3 - which is precisely what we have used in our 
formulae in the previous section. Specifically, while the Coulomb energy does not 
depend on the twist direction, the non-abelian energy is minimized when the largest 
(x 2 ') of the instantons is oriented in the lowest-potential direction X3. According to 
Eqs. ED and (l92l) . the instantons are larger in the direction of the twist than in 
other directions, hence the non-abelian energy 


E na = N C XM 1 + a 2 Ad ' 2 + 


- a 2 M 2 
2 


+ 


a 4 </>(2-7r — 4>) 
AD 5 


(M' 2 (n\ +n 2 ) + M 2 ri|) + 


0(a 4 M 4 ; 


(97) 


is minimized for n = (0,0, ±1), equivalently, 517(2) twist in the 73 direction. 

Also, Ad' > Ad reduces the equilibrium size of standalone or far-apart instantons 
from (1551) to 


, = (1/5 ) 1 / 4 

'° A 2 M 2 (M ' 2 + \ M 2 ) 1 / 4 


7 (for Ad' » M). 
VA MM' 


(98) 


Consequently the anisotropy, apart from stabilizing the chain in the X\. X2 direc¬ 
tions, provides an additional control over the equilibrium size of instantons. This 
will allow for an analytical analysis of the phase transitions below. 

Before we start the stability analysis we must describe the expected instability 
mode in the proper ADHM terms, i.e. in terms of T^ n matrices and Y£ vectors. 
Moving the instantons’ centers in 24 , 2,3 directions without changing their X 4 lo¬ 
cations along the chain, or any radii, or 517(2) orientations, as compared to the 
ansatz CD means keeping 


Y£ = (0,0, asin(n0/2), acos(n0/2)) 


r mn — Dn X S mn , 


(99) 
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exactly as for the straight chain but changing the rj^J ’ 2,3 matrices 
ran ^ ran [straight chain] + £r^ n , r^ n 

^ mn 5 ^ mn ^mn > 

( 100 ) 

in a manner that preserves the self-duality equations (l67l) . The solution correspond¬ 
ing to the displacement in the Xj, direction reads 

ST^l = 0, 6r 3 mn = S mn xSX 3 [n], (101) 

and the leading instability is expected to be the zigzag (figure [5]): 

5X 3 [n} = e x (-1)". (102) 


After defining the ADHM matrices we make a Fourier (Nahm) transform to 
map the infinite-dimensional matrices to differential operators on a circle. In the 
following it will be natural to combine wave functions ip{0) and i/j(6 ± tt) on the 
circle into a two-component wave function. Two-component functions provide a 
natural description for a two-layered chain. For some particular cases of the twist 
angle <f> = 0 or tt all expressions take particularly compact form. More generally 
for n layers n-component functions should be used instead. Specifically, we choose 
two-component functions with the following boundary conditions 


*(*) 


/ + ij){6 + 7r)\ 

W ( 0 ) - V >(0 + tt )/ ’ 


*(tt/ 2) = E 3 <F(-7r/2), (103) 


where Ei, E 2 , E 3 are Pauli matrices acting on the two components. In the space of 
such two-component wave functions, T 3 becomes the operator 


r> =exE 3 + ^e(-|<fl<|)xE 1 + ^(i-e(-|<9<f))xl, (104) 


where 0 is the step-function, i.e. 0 = 1 if —</>/2 < 9 < <f>/2 and 0 = 0 otherwise; 
and we restrict cj> to the values 0 < (j) < n. For cf> = ir the latter expression takes the 
form 

2 

T 3 = x Sr + e x S 3 • (105) 

Matrix T = <g) Y^ in the two-component notation becomes 


T(9) = ^(l + Y l )(s(9-^)+6(9 + ^)y 0 < 0 < tt . (106) 

For <f> = tt the latter can just be written as T = na 2 5(9±n/2) x 1, cf. equation (|83]l . 

Next we plug the data above in the formula for the operator L CD- To calcu¬ 
late the determinant of this operator, one needs to calculate its Green’s function 
L~ 1 {6 1 9q) and consequently the traces (jHTjl . As in the example of the straight chain 
the naive determinant is divergent and needs to be regularized, which is done by 
extracting an infinite but constant prefactor. One can reconstruct an analytic ex¬ 
pression for det(T) with any value of (j). However only for <j> = 0 and tt the expression 
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is compact enough to present in the paper. Accidentally </> = n is the minimal en¬ 
ergy configuration of instantons for both the straight chain at low density and in the 
zigzag phase provided that the zigzag amplitude is not too largef) In the remainder 
of this part we discuss the case <fi = n. One finds that the regularized determinant 
reads 


det(L) 

const 


( cosh TT 

2 

7r a 

sinh ?£) 

fcosh^ 

na 2 

sinh — 

\ D 

+ Dr± 

D ) 

V D 

+ Dr 2 
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r\ — (no 1 ID) 2 . nr i . nr 2 
— smh —— smh - 
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2rir 2 

nr x nr 2 

cosh —— cosh —— 
D D 


D 


D 


— cos ■ 


+ sin v cos 


— sin v cos ■ 


— sin v. 


7TCC4 

TT 

nxi 

IT 


2 cosh 


2 cosh 


nr i 

TT 

nr 2 

TT 


2ri r 2 

,2 


r 2 sinll ^l sinh 


D 


na . nri 

—— smli ——- 
Dr i D 

na 2 . nr 2 

—— smh ——- 

Dr 2 D 


2nx± 

TT 

nr 2 

TT 


(107) 


where 


r l ,2 = x l + x l + (x 3 =F 


b e = T 4e 2 + (na 2 /D) 2 , 

2e . 2e 

v = arctan ——— = arcsm — . 

7ra z /D b e 


(108) 

(109) 

( 110 ) 


Naively, we would expect the zigzag deformation to have no effect on the width 
of the instanton chain in the x\ and x 2 directions, while the width squared in the 
x 3 direction should increase by e 2 . And indeed, this is precisely what happens for 
any instanton radius and lattice spacing: evaluating (18911 . we obtain precisely 


2 2 2 4 

/ 2 \ / 2 \ a / 2 \ a n a 2 

(xi) = (x 2 ) = (xg) = — + + e . (Ill) 

Consequently, the non-abelian energy of the zigzag is, cf. (EHD, 

E na = N C XM {{\M 2 + M' 2 )a 2 + M 1 ^ + M 2 e 2 ^ 

= £ NA [e = 0] + N C XM 3 e 2 . (112) 


To calculate the integral in Eq. m and have workable expression for the 
Coulomb energy, we again resort to the approximation D a, where analytical 


f As we shall see in Sec. [5] this is true only under the assumption M' ^ M. Furthermore, it 
turns out if the zigzag geometry is relaxed there are other configurations favored over zigzag. See 
section 15.31 
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calculation is possible. For small, widely-separated instantons, r\ ~ a <C D, but 
finite e/D , we approximate 


det(L) « ( (a 2 + r\ + x\) + 


7T 


12D 2 

7ra;4 


(rf + 2a 2 r 2 - x\) + 0(a 6 /D 4 ) 

,2 


X I cos ■ D 

7rr 2 . . 7Tr 2 

-x - sinn —— 

8D 2 e 2 D D 


cosh 


7rr 2 

“d” 


7ra . 7rr2 

-— sinh —— 

2r 2 D D 


7r 2 a 4 


and consequently obtain 


E C = 


Nc■_ 
AM 


1 

5a 2 


3t r 2 
80Z3 2 


+ 


7r 2 tanh(7re/D) 


80D 2 


e/D 


+ 0(a 4 /D 6 ) 


(113) 

(114) 


Remarkably, the e-dependent part of this Coulomb energy is precisely five times 
smaller than the energy of point charges in a similar zigzag formation (1581) . Mathe¬ 
matically, this fivefold reduction stems from the last term in equation (11 1311 . which 
accounts for the interference between the instantons. It is not clear why the in¬ 
terference between well-separated small instantons has such a drastic effect on the 
Coulomb energy of the zigzag. Anyhow, the net energy cost of a small zigzag defor¬ 
mation e<dis 


Aifnet 


AFJna + AAq 


N C XM 3 e 2 


N c 

AM 


240.D 4 


600.D 6 


0(e 4 /D 6 ) 


(115) 


This cost is positive 
large lattice spacings 


and hence the straight chain is stable 


D > D c 


7 T 1 

^240 My/X ' 


for all sufficiently 


(116) 


For smaller lattice spacings D < D c , the energy function (1 115 D has a negative coeffi¬ 
cient of e 2 but positive coefficient of e 4 . Thus, for D < D c the straight chain becomes 
unstable and there is a second-order phase transition to the zigzag configuration. 
For lattice spacings just below critical, the zigzag parameter e is 


(e) m ±— y/D c (D c - D). (117) 

7T 

For smaller lattice spacing (but larger than the instanton size), the zigzag parameter 
satisfies the same equation as in the point charge limit (l60l) with the new D c . 

In units of the equilibrium instanton size (f98l) . the critical lattice spacing (11161) 
is 


^crit 


' M r 


\M 2 


a « a n 


v^48 V 


M 2 


1/4 


M’ 

, for M'>M. (118) 

M v 1 


For a highly anisotropic gauge coupling with M' M, the critical spacing is 
much larger than the instantons, which justifies our approximations. However, for 
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near-isotropic couplings with M' « M, the critical spacing is only -D cr it ~ l-32ao, 
and in this regime we cannot be sure our analysis of the zigzag instability is even 
qualitatively correct. That is, we are not sure that there is a transition from a 
straight instanton chain to a zigzag for M' ss M , never mind any details of such 
transitions. 

In this section we have presented analytical solutions for two types of instanton 
lattices: the straight chain and the zigzag with anti-ferromagnetic orientation of the 
instantons. More specifically, we have spelled out the result for the determinant of 
the operator L, which can be readily used to evaluate the instanton density via 
Eq. m- The result allows to compute easily the non-abelian contribution to the 
energy of the configurations, however the integrals giving the abelian Coulomb part 
already require some approximations, such as the small instanton limit a <C D. 

The configurations considered here seem to approach the limits of what can be 
done exactly. More complicated instanton lattice configurations require simplifying 
assumptions already at the level of the ADHM ansatz. This will be done in the 
next section, where the two-body force approximation will be implemented. The 
approximations allow to find many new phases, including the ones with non-abelian 
twist angles. To see them however, one needs to further increase the anisotropy of 
the system. In all cases, however, we will assume the small instanton limit. 

5. Approximations and generalizations 

In this section we will review other interesting phases of the holographic nuclear 
matter by exploring a shortcut around steps (1), (2), and (3) in section lT3l Following 
the idea of Ref. EDI we will argue that when the distances between the instantons 
are much larger than their radii, the interactions between A baryons are dominated 
by the two-body forces, 

» \ Y, MX&-XZ, m,n) (119) 

m,n=l,...,A 

m^n 

for a manageably simple function F 2 of the two instanton’s positions X! r J n n and 
orientations m and n. Consequently, we can minimize the net interaction energy 
over the entire moduli space of the multi-instanton system using a simple numerical 
simulation: Starting from a random set of instanton positions and orientations, use 
the steepest descent algorithm to find the nearest local minimum of the net energy; 
repeat this procedure for different random starting points to find other local minima; 
eventually, find all the local minima, compare their energies, and identify the global 
minimum. 

In the numerical simulation in Ref. EDI the lattice geometry for the instanton 
positions was assumed, but the numerical simulation was used to find the lowest- 
energy pattern of their orientations. Once all the patterns are known, they can be 
used as ansatz’s (with a few parameters) for which the net energy can be calcu¬ 
lated as analytic functions of the parameters. In the boundaries between different 
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orientation patterns are mapped much more accurately than just in the numerical 
approach. In this section we will review the main results of Ref. 1301 

5.1. Two-body forces 

In real-life nuclear physics, besides the two-body nuclear forces due to meson ex¬ 
changes, there are significant three-body forces, and presumably also four-body 
forces, etc., 

A i 

^nucleus = ^ff lbody (n) + - £ H 2hody (m,n) 

n= 1 different 

m,n=l,...A 

+ l H 3hody (£,m,n) + ••• (120) 

different 

£,m,n=l,...A 

where n stands for the quantum numbers of the n — th nucleon. Likewise, in the 
holographic nuclear physics interactions between multiple baryons include two-body 
forces and also three-body, four-body forces etc. Even in the classical infinite-mass 
limit, where the holographic baryons become static instantons of the SU(Nf) gauge 
fields in 4 + 1 dimensions, the potential energy (due to five-dimensional curvature 
and due to Chern-Simons interactions) of an A-instanton system has form 

AE(l,2,...,A) = £ total - AXN C M 

a ^ i _ i 

= ^ £ (1) (n) + - ^2 E {2 \m,n)+- ^ £ (3) (l, to, n) + • • • (121) 

n= 1 different different 

m,n=l,...A £,m,n=l,...A 

with significant three-body, etc. terms. 

What about the relative magnitudes of the two-body, three-body, four-body, ..., 
interaction terms? When the baryons are tightly packed so that their instanton cores 
overlap and merge, we expect all the n-body forces to have comparable strengths. 
But in the opposite low-density regime of baryons separated by distances much 
larger then their radii, the two-body forces dominate the interactions, while the 
multi-body forces are smaller by powers of (radius/distance) 2 . In this subsection 
we briefly summarize the proof of this observation given in Ref. 1301 

We review the simplest case oiNf = 2, although the proof can be extended to the 
general SU(Nf) case as well.® First let us consider intermediate-range distances 
between the baryons, 

( a ~M7x) <<: \ Xm ~ Xnl « ~M ’ (122) 

which allows us to treat the five-dimensional holographic space as approximately 
flat. The holographic baryons are instantons of the SU( 2) magnetic fields, which 
source the 1/(1) electric and scalar fields via CS and tr (^F^F^) couplings. When 
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those instantons are small and separated from each other by large distances \X m — 
X n | ~ D a, their interactions come from two sources: 


(1) Direct Coulomb repulsion (electric—scalar) between nearly-point-like abelian 
charges in 4 + 1 dimensions, 

N P . X ^ 1 


A E' 


direct 


E 


4AM ^ \X^-X^\ 2 ’ 

m^tn 1 1 


(123) 


which is a manifestly two-body interaction. 

(2) The interference between the instantons. The latter changes the distribution of 
the instanton number density in space, 




I{ x ) — 227J-2 ( a ')^per (* E )) ) 


(124) 


IA instantons ^ ^ ji 


standalone 


(®) + AT 


interference 


(*), (125) 


which in turn changes the self-interaction energy (Coulomb and non-abelian) 
of each instanton by an amount comparable to (112311 . 


Within an instanton, i. e. at 0(a ) distance from some instanton’s center X%, the 
interference from the other instantons should be relatively weak, 

,,2 

A ^interference „ x /standalone/^ (126) 

D 1 

so there should be some kind of a perturbation theory for it. At the first order of 
such perturbation theory, the AJ arises from interference between the un-perturbed 
standalone-like instantons, so we expect it to be a sum of pair-wise interferences 
from the other instantons, 

AJ l S t order (;cnearXji) = /standalone^) ^ , (127) 

m^n 

where 

x£l( x ) ~ °(\x n -x m \*) (128) 

and depends only on the instantons and i.e. on their positions, radii, and 
orientations, but not on any other instantons. 

At the second-order, we expect to include the interference between the first-order 
A I and the additional instantons, so at this order we obtain three-body effects, 

A / 2nd order (a . nearAn ) = /sitandalone^) ^ ^J X ) (129) 
















40 V. Kaplunovsky, D. Melnikov, J. Sonnenschein 


Likewise, the higher orders may involve more and more instantons, but the magni¬ 
tudes of such high-order interference effects are suppressed by the higher powers of 
(a 2 / D 2 ). 

Now consider the Coulomb self-interaction of the instanton #n, 


E6 elf (n) = 


N r . 1 

X 


4AM pi 

where p n is the instanton’s effective charge radius, 


1 


x\ ,X 2 near X n 


,4 , 4 I(x 1 )I(x 2 ) 

CTX 1 <TX 2 -ry 

Fi - X 2 J 


A standalone instanton has 


1 4/5 

Z2 - XT 


(131) 


(132) 


(133) 


Pn K 

but interference from the other instantons should change this radius by small 
amount of similar relative magnitude to A///® tandalone , thus 

1 4/5 1 

Pi a 2 n + n ’ " ~ D 2 ’ 
which changes the instanton’s Coulomb self-interaction energy by 


(134) 


A E^\n) = 


N r . 


X Nn 


N r . 


(135) 


4AM '' AMD 2 

Note that this effect has a similar magnitude to the direct Coulomb repulsion (I123D 
between the instances. 

Moreover, the charge radius correction A n is linear in the A/ lnterference at x near 
the X n , hence in light of Eq. (11271) . the leading-order contribution to the A„ is a 
sum of pair-wise interferences from the other instantons, thus 

A™ = ^ AW m + 0(a 2 /D 4 ), (136) 

m^n 

where each Adepends only on the instantons and #m. Consequently, the 
leading effect of the interference on the net Coulomb self-energy of all the instantons 
has form 

A int £^ elf = D£ elf [interfering] - Eg* 1 [standalone] 

^ y A w " 2 

/ j n,m 


4AM 


O 


( N c a“ 

I -x 


n,m=l,...,A 
n^m 


{AM D 4 


(137) 


where the leading terms act as two-body interactions between the instantons. 

Besides the Coulomb energy, the non-abelian energy is also affected by the in¬ 
terference between the instantons, 


A mt D NA = N c AM i / cfx x 4 x A/ 


interference 


(a), 


(138) 
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but this time the integral should be taken over the whole four-dimensional space, 
including both the instantons and the inter-instanton space. Indeed, one can show 

thaP! 


at 0{a ) distance from an X n , A I 
at 0(D) distances from all the X m , 


1 

jitd? > 


(139) 


A I 


D s ’ 


and both kinds of places make 0(a 4 /D 2 ) contributions to the integral (|138[) . More¬ 
over, in both kinds of places, the leading terms in the A I(x) is a sum of independent 
two-instanton interference terms, 

^interference^) = * ^ l^x) + Subleading. (140) 

n^m 


Although our heuristic argument (based on A/ mterference <C / standalone ) for such 
decomposition near instanton centers does not work in the inter-instanton space, 
we shall show later in this section that the decomposition ( 11401 ) works everywhere 
in the four-dimensional space. Therefore, the non-abelian interactions between the 
instantons due to interference are also dominated by the two-body terms. 

Let us now relax the un-necessary assumption of intermediate-range distances 
between the baryons and consider what happens at longer distances D ~ 1/M. In 
this regime, the curvature of the fourth space dimension — and especially the z— 
dependence of the five-dimensional gauge coupling ( 1351 ) can no longer be treated 
as a perturbation. Consequently, the magnetic fields of an instanton or a multi- 
instanton system are no longer self-dual in the inter-instanton space. Therefore, 
the interference between very distant instantons is no longer governed by the self¬ 
dual ADHM solutions. Instead, we must work it out the hard way: figure out how 
the magnetic fields of an instanton propagate through the curved 4+1 dimensions 
towards the other instantons, and then find out how such fields disturb the other 
instantons’ cores. 

Fortunately, we do not need a hard calculation to see that at large distances from 
an instanton its magnetic fields are very weak. Indeed, even in flat five-dimensional 
space the fields weaken with distance as A“ ~ a 2 /r 3 (in the IR-safe singular gauge), 
so at r>a they are so weak that the field equations become effectively linear. In the 
curved space, we may decompose the weak hve-dimesional gauge fields into four¬ 
dimensional mesonic fields, hence at large distances r > 1/M from an instanton, its 
fields become Yukawa-like 

2 2 
A;(r,z) ~ ^£tf*(;0**Cw)e-"» r < (141) 


Consequently, at the location X!/ h of any particular instanton, the background fields 
from the other instantons are very weak, and their effect on the instanton //n itself 
can be adequately accounted by the first-order perturbation theory. In other words, 















42 V. Kaplunovsky, D. Melnikov, J. Sonnenschein 


the effects of other instantons ^ n on the instanton are weak and add up 
linearly. For the self-interaction energy of the n th instanton, this means 

E(n) ss .E(standalone) + Ea m E(n), (142) 

m^n 

where the second term gives rise to two-body interaction energies 

^interference ( n > m) = A m E(n) + A n E(m). (143) 

We expect the two-body terms to be rather small: in addition to the usual 
\/\X m — X n \ 2 factors from five dimensions they should carry Yukawa-like exponen¬ 
tials e~ mr (or rather sums of such exponentials), but these are the leading inter¬ 
actions due to interference. The three-body or multi-body interactions follow from 
higher-order perturbation by very weak fields (11411) , so they are much smaller than 
the two-body interactions. 

As to the direct Coulomb interactions between the holographic baryons via elec¬ 
tric or scalar fields, at large distances, \X m — X n \ ~ 1/M, they also decompose into 
sums of Yukawa forces. Moreover, since the scalar mesons generally have different 
masses from the vector mesons, the attractive potential due to scalars may have 
a different r dependence from the repulsive potential due to vectors. Thus, for a 
right model, the net two-body force between two holographic baryons may become 
attractive at large enough distances between them. But regardless of the model, 
the direct Coulomb interactions are always manifestly two-body for Nf = 2, while 
for Nf > 3 the multi-body terms exist but become very small at large distances 
between the baryons. 

Thus, 

^direct = ^ E overlap (m,n) x F (2) (|A m - X n \), (144) 

mpn 

where the precise form of the potential V^ 2 \r) is model-dependent, but the two- 
body form of the direct interactions is quite universal. 

To conclude this part, let us demonstrate how the multi-body expansion works 
by evaluating the non-abelian interaction energy of an A-instanton system. Due to 
the nice tricks with integration by parts of the instanton density in Eqs. (l87l) (1M1) 
the expansion can be worked out quite easily. This is not the case with the Coulomb 
energy and we refer the reader to Ref. [301 for the full details of the proof. 

We would like to implement the large distance expansion for the solution of the 
ADHM equations for the off-diagonal elements T^ = al^ n , see section l4~Tl Indeed 
for large distances \X m — X n \ ~ D between the instantons, D a„, we may solve 
the ADHM equations, together with constraints (171)1) as a power series in a 2 /D 2 : 


= a (1) + a (2) + a (3) 

mn ^fimn 1 u /imn 1 imn 


( 145 ) 
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where 


a 


(i) 

fimn 
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( 2 ) 

gmn 


a ( 3 ) 

p,mn 


(X m 

\x m - 


-X„) 




X n | 2 


x i a m a n tY{yl n y n (-iT a )) 


(X m X n ) v 
\Xm - X n \ 2 


X V n a rr (1) rv (1) = 

£^m,n 


= 0(a 2 /D ), 
0(a 4 /D 3 ), 


^1% ( X m -X n ) v 

\x m - x n \ 2 


X V n a rv (1) rv (2) 
£^m,n 


0(a 6 /D 5 ), 


(146) 

(147) 

(148) 


and so on. Note that for each off-diagonal matrix element, the leading term aj,™ 
in this expansion depends only on the instantons jfcm and Jfn (that is on the 
positions, radii, and orientations of only these two instantons), while the subleading 
terms afin, ■ • ■ involve additional instantons. 

Let us plug this expansion in the expression for the second moment of the 
instanton number density, which computes the correction to the non-abelian energy 

4 

A£ na = X>4a, (149) 

p=l 

where the sum goes over the contributions of all dimensions and 
A E^l = \N C MMl J d 4 xl(x) x (x p ) 2 = XN c MM 2 (tr(r»T») + t tr(T)) 

= A N C MM 2 (j2{( T nn) 2 + \Tnn) + E( r --) 2 

\ i=n m^n 

= XN c MM 2 J2((XZ) 2 + \a 2 n ) + A N C MM 2 £ «„) 2 . (150) 

n m^n 


The above correction appears due to curvature M 2 in the direction cf. Eq. m 
for the five-dimensional coupling. Obviously, the first sum on the last line here is 
the sum of individual instantons’ potential energies due to their radii and locations 
(relative to the x M = 0 hyperplane), while the second sum comprises the interactions 
between the instantons. 

Irrespectively of the details of the instanton configuration, the small instan¬ 
ton expansion (11451) tells us that the interaction term in Eq. (11501) will contain 
only the two-body terms in the leading order in a 2 /D 2 . The leading order term is 
0(XN c M 3 a 4 / D 2 ), while the multi-body term will be 0(XN c M 3 a 6 / D 4 ), or weaker. 
Specifically, 


A-^na (2-body; to, n) 


X N c MM 2 a 2 m a 2 n 
2\X m -X n \ 2 




(151) 


where N mn is the 3-vector, which contains all but y-th component of the unit 4- 
vector 


N " = ( N 1 N 2 N 3 N 4 ) = 

mn V mn ’ mn ’ mn ’ mn J 


vu _ vv 

- / *-m 

\x n — x m \ 


( 152 ) 
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Let us also quote the result for the Coulomb energy obtained through a tour 
de force calculation in Ref. 1301 In the multi-body expansion of the net Coulomb 
energy, 

£totai = + l J2 E c h ° dV ( m ’ n ) + , (153) 

n m,y£n 


the one-body and two-body terms are found to be 


£^ body (n) 


E* hody 


(m, n) 


5AM a 2 ’ 

No 1 ( 1 + l(“k 

2AM \X m — X n \ 2 \ 5 \ a 2 


(154) 

4 ) (tr 2 (yly n )~2)), 

a m J J 

(155) 


For small instantons distant from each other, the one-body potential energies for 
the instanton radii a,; are much larger than the two-body etc. interactions between 
different instantons. Consequently, in the minimal-energy or near-minimal-energy 
configuration of the multi-instanton system, the instanton radii will be close to the 
equilibrium radius of a stand-alone instanton, 

a n = a 0 + 0(a 3 /D 2 )z , (156) 


In particular, for the most general case here of non-zero M 2 , M 3 and M 4 , cf. Eq. ([55]) , 


a 


4 

0 


2/5 

A 2 M 2 (M| + Mf + Mf) ’ 


(157) 


Plugging this equilibrium radius into the combined non-abelian and Coulomb two- 
body interaction energy one obtains 


NE^ t {m,n) = 



where 


C, = 


M 2 

±v± /l 


M 2 


M 2 


Mf ’ 


Ci + C 3 + C 2 = 1. 


(158) 


(159) 


For the original holographic setting with Mi ^.3 = 0, while M 4 = M, this formula 
simplifies to the Kim and Zahed formula!^ 


E 2hody (m,n) 


2 N c 1 
5AM \X m - X n \ 2 





(160) 

Note that the expression inside the square brackets is always positive, so the two- 
body forces between the instantons are always repulsive, regardless of the instantons’ 
SU (2) orientations. However, the orientations do affect the strength of the repulsion: 
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two instantons with similar orientations repel each other nine times stronger then 
the instantons at the same distance from each other but whose orientations differ 
by a 180° rotation (in SO( 3) terms) around a suitable axis. This fact will be at the 
core of our analysis of instanton crystals in the remainder of this section. 


5 . 2 . Phases and transitions 

The perturbative analysis of instanton solutions and their interaction energies, based 
on the multi-body expansion discussed above, allows to vastly expand the search 
of available phases of one-dimensional lattices (both straight chains and zigzags). 
In particular, numerical algorithms can be implemented for a search of the ground 
state for a given set of parameters. Indeed we do not know a priori which choice 
of instantons’ moduli parameters, such as positions and orientations, will minimize 
the net energy. Actual realization of the exact solutions presented in section 0] so far 
can only be conjectured. In view of a large number of possibilities the right strategy 
seems to be to start from a random values of the moduli and let them evolve and 
relax numerically to a local energy minimum. Comparing the energies of the local 
minima a true ground state can be selected. 

The strategy outlined here was partially realized in Ref. 1301 where the positions 
of instantons was assumed to fill a one-dimensional chain, or a zigzag, but the 
orientations were allowed to take any values in the numerical simulations. The 
numerical relaxation have indicated a set of minimum energy configurations, from 
which the true ground states were selected and preliminary phase diagrams were 
obtained. Given the ground states, the analytical ansatzes were constructed for a 
more accurate detection of the location of the phase boundaries. Let us review the 
main results of that investigation. 

Following Ref. [30], we assume that the holographic direction z = x 4 , and there is 
an additional anisotropic curvature in the directions x 2 and a ; 3 so that the instantons 
are stabilized in a one-dimensional array along the direction x 1 . This corresponds to 
the choice of five-dimensional gauge coupling as in Eq. (l52l) with M 4 = M. We will 
also assume M 2 <C M 4 , while M 3 will be allowed to take any values M 2 < M 3 < M 4 . 
In such a setup the instantons will line-up along x 1 at low densities and the leading 
instability will be the zigzag deformation in the ±x 2 direction. We will also see that 
the lowest energy configuration will be sensitive to the anisotropy ration M 3 /M 4 . 

First, let us review various possibilities for the straight chain, 

X% = (nD, 0,0,0), neZ, (161) 

where we will assume that the lattice spacing D is much larger than the instanton 
radius a. The curvature creates the following effective one-body potential for the 
instanton centers X n , 

A E {1 \n) = A N C M(M^{X 4 ) 2 + M 2 (A 3 ) 2 + M 2 (X 2 ) 2 ) + 0(N c XM 5 X 4 ). 

(162) 
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First, consider the regime M 2 < M 3 <C M 4 . In this limit, the M 2 and M 2 
parameters give rise to the x 2 and £§ terms in the one-body instanton poten¬ 
tial (11621) , but their effect on the instanton radius a or the two-body forces between 
the instantons may be neglected (compared to the effect of the M|x| term). There¬ 
fore, the net two-body forces between the instantons remain approximately as in 
Eq. (11601) . which for the one-dimensional lattice geometry (11611) of the instanton 
centers \X m - X n \ 2 = D 2 {m - n) 2 and for N mn = (N^ n , N^ n , N^ n ) = (±1,0,0) 
may be summarized as 


A E 


( 2 ) = 
int 


N c y^ 
5AM D 2 ^ 

m^n 



1 

2 


+ tr 2 {ylny n ) + tr 2 (yl,y n (-iTi)) 


(163) 


As we noted in the end of section l5Tl to minimize this energy, each pair of 
instantons m and n wants to have the relative orientation being a rotation through 
a 180° angle, in 50(3) terms, around some axis perpendicular to x 1 . In other words 
it wants y' m y n to be a linear combination of zr 2 and zr 3 . However, this cannot be 
achieved for all instanton pairs at once: minimizing the energies of the (n, m) and the 
(n,£) pairs will not produce a minimum of (m,£). Instead, we will first minimize 
the energies of the nearest-neighbors and then look what happens with the less 
expensive pairs of instantons. 

Thus, we want 


Vto : ylnVm +1 = cos ^ m (zr 3 ) + sin (^ 2 ) for some angle ip m , (164) 


which has the following most general solution (modulo a common SU(2) symmetry): 


y n = exp(i^„Ti) (*T 3 ) n 


±[cos</>„ 1 + sin <j> n (zn)] for even n, 
Tfcos^n (zt 3 ) + sin cj) n (zt2)] for odd n, 


(165) 


for some angles (f> n . In particular, = (—l ) m — 4> m ). 

The big family of solutions (I165D turn out to be completely degenerate! Indeed, 
for all sets of (j > n , all instanton pairs (m, n) with odd m — n have minimal energies, 
while pairs with even m — n have maximal energies. Consequently, regardless of the 
angles 4 > m , the net energy (11631) per instanton of the chain is 


Flint _ X c j 1 1 9 1 

(1-inst) “ 5A MD 2 1 2 ^ l 2 + 2 ^ l 2 

^ Y odd^ even £7^0 

In a generic lowest-energy configuration (I165D . the instanton orientations span 
an 50(2) x Z 2 subgroup of the 50(3) = 51/(2) x Z 2 of the y n , which corresponds 
to the rotational symmetries of a cylinder — rotations through arbitrary angles 
around the aq axis, and 180° rotations around axes perpendicular to aq. The y n 
alternate between the two types of rotations, but apart from that they generically 
do not follow any regular patterns. 

However, the family (11651) also contains some regular patterns in which the 
orientations y n (modulo sign) follow a repeating cycle of finite length p; moreover, 


IT 2 N r 


10A MD 2 


(166) 
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the values of y n span a discrete subgroup of the cylindrical symmetry SO( 2 ) x Z 2 . 
Here are some examples: 


• The anti-ferromagnetic chain , with two alternating instanton orientations: 


2/evenn — d=l, 2/odd n — T 2 T 3 . 


(167) 


In this configuration, for cj> n = 0, the y n (modulo sign) span a Z 2 subgroup of 
the 50(2) x Z 2 . 

• Period-four configuration spanning the Klein group of 180° rotations around 
the three Cartesian axes@ 


(168) 


Vn= 0 (mod 4) — dll, 2/n= 1 (mod 4) — dzT 3 , 

Vn= 2 (mod 4) dlTl, Vn=3 (mod 4) dz72. 

• Period-2fc, 2k = 6 , 8 ,10,... configurations spanning prismatic groups r L\ z x Z 2 : 


7T71 . 7m 

2/even n — COS ~7y^ L H T Sin TyT (i T\ ), 

n(n — 1 ) .. . 7 T (n — 1 ) 

2/odd n = COS - — - (*r 3 ) + sin -—- (it 2 ). 


(169) 


2k v ' 2k 

• Period-2fc, 2k = 6,8,10,... configurations spanning dihedral groups D 2 k, for 


(j> n = n x (7r/2fc), ie. 


2 /e 


7 rn 7rn 

= C0S 2 A 1 + Sm 2 fc ( * Tl) ’ 


7TTI 7TTI 

2/odd n = cos —(zr 3 ) + sin —(ir 2 ). 


(170) 


There is a wider class of regular configurations, we shall call them link-periodic, 
in which the y n themselves are not periodic, but the relative rotations 2 /i 2 /n+i be¬ 
tween nearest neighbors follow a periodic pattern. In terms of the <f> n angles, this 
corresponds to periodic differences ^ n +i — <p ni for example 


mp 



ip + 9 for even n, 
ip — 6 for even n. 


(171) 


For rational <p/n this pattern produces a periodic array of instantons’ orientations 
y n , e.g. the dihedral cycle for ip = n/2k and 9 = 0 , or the prismatic cycle for 
ip = 7 t/2 k and 9 = —p/2. For irrational p's the orientations y n themselves do 
not have a finite period; instead, they wind irrationally around the cylinder group 
5/7(2) x Z 2 . 

The large degeneracy of the one-dimensional configurations of instantons can 
be lifted if the condition M 2 <C M 4 is relaxed. For M 2 , M 3 ~ M 4 , there is a unique 
lowest-energy straight chain configuration, namely the link-periodic array (11711) . 
whose p and 9 parameters depend on the M 2 /M 3 ratio. To see this one can use 


g As a subgroup of 50(3), the Klein group is abelian and isomorphic to Z 2 x Z 2 . But its cov¬ 
ering group in SU( 2) is non-abelian group and isomorphic to the group of unit quaternions 
=tl, =b i, d zj, d zk. 
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Eqs. (11501) and (11581) for the one-body and two-body interaction energies in sec¬ 
tion [5T] 

Due to a larger number of parameters it is less trivial to minimize the interaction 
energy. In particular, it is not intuitively obvious how to balance the energy cost 
of nearest-neighbor interactions versus next-to-nearest neighbors and more distant 
instanton pairs. To find the ground state, in Ref. [30] a numerical experiment was 
performed using a lattice of 200 SU( 2) matrices y n . In each run, y n ' s were assigned 
some initial random values, and then they were alowed to evolve towards a minimum 
of the energy function (|158D via the relaxation method. 

The bottom line of the numerical simulations in Ref. (30] is that the ground 
state of the one-dimensional instanton lattice always has a link-periodic instanton 
orientations (11711) (up to a global symmetry, if any) for a periodicity angle <p that 
depends on the ratio of parameters M 2 /M 3 and not on the M 4 (as long as M 4 > 
AT) > M2). 

Once the ground state is known the exact dependence of the angle ip on the 
parameters M 2 and M 3 can be calculated analytically. Using Eq. (11711) for the 
ulnUn and the expression for the energy (|158l) one finds 

(St) = (t) 1 + & + ft - (ft - ft)cos(2ft)) 

- 2 ttM (c 3 + C 2 - (C 3 - C 2 ) cos (2 0)) + V (C 3 + C 2 )) ■ (172) 

Minimizing this expression with respect to the ip and 9 produces four degenerate 
minima, namely 


(1) 


7 r 


Cb 

<p = 

= + 2 

X 

c 2 + c 3 ’ 

(2) 


7 r 


Co 

<p = 

~2 

X 

c 2 + c 3 ’ 

(3) 

p = 

7 r 

~2 

X 

C 3 

c 2 + c 3 ’ 

(4) 

p = 

7 r 

= + 2 

X 

C 3 

c 2 + c 3 ’ 


however the last two minima are physically equivalent to the first two. In agree¬ 
ment with the numerical experiment, the one-dimensional instanton lattice has two 
degenerate ground states related by the p> —> —ip symmetry. The value of \ip\ is also 
independent from M 4 . 

Two particularly interesting M|/M| ratios need special handling, Mf = Mf and 
M 2 -C Mf. For M 2 = M 3 the background has a rotational symmetry in the x 2 ' 3 
plane. For the instanton chain, this translates into the 17(1) symmetry between the 
ir 2 and it 3 directions in the 517(2). Consequently, instead of two discrete ground 
states of Eq. (11731) we have a continuous family: 

7r 


9 = anything. 


(174) 
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For all configurations in this family, the instanton orientations y n (modulo signs) 
repeat with period four, while spanning the Klein groups Z 2 x Z 2 of 180° rotations 
around three mutually perpendicular axes. For 0 = — f the y n are spelled out in 
Eq. (11681) : for other values of 8 we have similar cycles related by the U( 1) C SU(2) 
symmetry. 

Finally, for the very asymmetric background with Mlj <C M|, the two ground 
states of Eq. (11731) become indistinguishable as <p —\ ±0. In this limit, the instanton 
orientations form the anti-ferromagnetic order (11671) . This is indeed the case we 
have considered in section 14.11 if we assume M 3 = M 4 = M' and M 2 = M -C M'. 

The analysis of the zigzag configuration proceeds similarly. We remind that the 
original chain was in the X\ direction, while by virtue of the condition M 2 < M 3 < 
M 4 the most likely first step in the transition between one-dimensional and two- 
dimensional instanton lattices is the zigzag in the ±22 direction (Fig. [5]). The ADHM 
ansatz for the zigzag appeared in Eqs. (l99ll . (11011) and (11021) . It is parameterized by 
two parameters D and e; we shall refer to D as the lattice spacing and to e as the 
zigzag amplitude. 

In the analysis of the zigzag transition we shall use the small instanton approx¬ 
imation, a <C D. It was demonstrated in section l4~2l that for consistency one need 
to impose the condition M 2 <C M 4 z = M. Indeed in such a case the critical spacing 
for the zigzag transition, will respect the a <C D approximation, as in Eq. (11161) . 
The parameter M 3 , instead, is allowed to take any values M 2 < M 3 < M 4 . Thus, 
for the search of minimum energy configurations, one can evaluate the net inter¬ 
action energy per instanton using Eq. (11581) . with negligible C 2 ~ 0. For a zigzag 
where all the instantons lie in the x 1,2 plane and hence all the N^ n have form 
Nm n = (*,*,0,0), Eq. (11581) can be rearranged as 

[zigzag] = N c y. Q z {m,n ) 

(1-inst) 5AM ^ \X m - X n \* ’ 1 j 

m^n 1 ' 

where 

Q z {m,n) = \ + tr 2 (yl n y n ) + C 3 ^ tr 2 (y] n y n {-iT a )) 

< 2 = 1,2 

+ (1 - 2C 3 ) tr 2 (ylnVni-iT ■ N mn )). (176) 

Note that in general, that the force between two instantons is not central — it de¬ 
pend not only on the distance between the instantons and their relative orientation 
ylnUn °f their isospins, but also on the direction N mn of their separation X m ^X n in 
space. However, in the background with M 3 = M 4 , which was considered in Ref. [29' 
(albeit in different notations) — C 3 = 1/2, the last term in Eq. (11761) goes away, 
and the force becomes central (but isospin-dependent). As a result, for M 3 = A/ 4 , 
the only minimum energy configurations found in Ref. 1291 corresponded to relative 
orientation of instantons generated by an abelian U( 1) in SU( 2). For M 3 ^ M 4 , 
as we shall see, the phase space is richer. Indeed, let us discuss the lowest-energy 
configuration of the orientations y n as a function of the e/D and M 3 /A /4 ratios. 































50 V. Kaplunovsky, D. Melnikov, J. Sonnenschein 


As in the case of the straight chain, it is hard to say a priori , which particular 
pattern, or patterns the relative orientations y/jjn -fi of the neighbor instantons 
will follow. It is thus useful to start with a numerical simulation, starting with 
completely random y n and letting them evolve towards a minimum of the energy 
function (11751) . In Ref. [301 such a simulation for different combinations of the e/D 
and M 3 /M 4 parameters was performed and the following four patterns of orientation 
was discovered: 


• the anti-ferromagnetic pattern (AF) of instanton orientations, in which the 
nearest neighbors always differ by a 180° rotation around the third axis, 


Vn 



for even n, 
for odd n, 


same y\y n+1 = ir 3 for all n ; (177) 


• another abelian pattern (AB), in which all nearest neighbors differ by the same 
17(1) C SU(2) rotation, but now the rotation angle is less than 180°, 

same yj l y n+1 = exp^4>T 3 ) for all n, 0 < <j> < 7 r; (178) 

• a non-abelian link-periodic pattern (NAl), in which the relative rotation be¬ 
tween nearest neighbors is always through a 180° angle, but the direction of 
rotation alternates between two different axes in the x 1,2 plane: one axis for 
the odd-numbered instantons and the other for the even-numbered. In SU( 2) 
terms, 


yl k y 2 k+i = ex P ( — n e ■ = ift e ■ t = +iAr 3 + *Br 2 , 

V \ ( 179 ) 

y\k+iV2k+2 = exp ( — n 0 ■ t j = in 0 -r = +iAri - iBr 2 , 
for some A,B ^ 0 ( A 2 + B 2 = 1). 

• another non-abelian link-periodic pattern (NA2). Again, the relative rotation 
between nearest neighbors is always through a 180° angle and the direction of 
rotation alternates between two different axes. However, this time the two axes 
no longer lie within the x 1,2 plane, thus 

yl k y2k+i = iAr i + iBr 2 + iCr 3 , 

(loll) 

V 2 k+iy2k+2 = i Ar i - i Br 2 - iCr 3 , 
where A,B,C^ 0 (A 2 + B 2 + C 2 = 1). 


Figure [G] shows the distributions of the four phases in the parameter space. 

Figure [ 6 ] also shows the order of the phase transitions between distinct phases. 
In particular, the transitions between the two non-abelian phases are second-order. 
(Notice, that the phase non-abelian NAl appears in two disconnected regions of the 
diagram, separated by NA2.) Likewise, the transition between the non-abelian NA2 
and the abelian anti ferromagnetic phase AF is also second-order. On the other 
hand, the transitions between the other abelian AB phase and the phases NAl and 
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M3/M4 


A 



Fig. 6. Phase diagram of the instanton orientation patterns in the zigzag parameter space. This 
diagram is obtained from the mostly-analytical calculation in Ref. 1301 The lines separating different 
phases indicate the order of the phase transition: a white line for a second-order transition and a 
black line for a first-order transition. 


NA2 are first-order. Finally, the triple point at (M 3 /M 4 ) = 1, ( e/D) ss 0.38 between 
the AF, NA2 and AB phases is critical. At that point, the transitions between all 
three phases are second-order. 

Thus far, we considered the zigzag amplitude e as an independent input param¬ 
eter. But physically e is a dynamical modulus whose value follows from minimizing 
the net energy of the multi-instanton system. In fact, all instanton center coordi¬ 
nates are dynamical moduli, which raises two problems. First, in some situations 
the lowest-energy configuration of the instanton lattice may be more complicated 
than a zigzag or a straight chain. Second, for some lattice spacings A, a uniform 
lattice of any kind — a straight chain, a zigzag, or anything else - - may be unstable 
against breaking into domains of different phases with different densities. 

By way of analogy, consider a fluid governed by an equation of state such as 
Van-der-Waals. Formally, this equation allows a uniform fluid to have any density 
(up to some maximum). But in reality, at sub-critical temperatures one may have 
a low-density gas or a high-density liquid, but there are no uniform fluids with 
intermediate densities. If we constrain the overall volume V of some amount of fluid 
such that its average density would fall into the intermediate range, we would not 
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get a uniform fluid; instead, part of the volume would be filled by the higher-density 
liquid while the other part by the lower-density gas. In the same way, if we fix the 
overall length L of the x 1 axis occupied by some large number A of instantons, 
their lowest-energy configuration is not necessarily a uniform lattice of some kind. 
Instead, for some average density p = A/L we we would have L split into domains 
of two different lattices of different densities. 

To keep the fluid uniform, one should control its pressure P rather than the 
volume V ; consequently, the preferred phase follows from minimizing the free en¬ 
thalpy G = E — ST + PV rather than the free energy F = E — ST. Likewise, for 
the one-dimensional lattice, we should control the net compression force F along 
the x 1 axis, rather than the net length L or the lattice spacing D. Also, we should 
minimize the free enthalpy of the lattice G = E — ST + LF, but since we work at 
zero temperature all we need is the ordinary enthalpy H = E + LF. Equivalently, 
we may minimize the non-relativistic chemical potential 

Gnon-rel E + LF F 

P = Prel - Mbaryon = - 1 - - A - = £ + ~ ■ (181) 

1 —0 J\. p 

(We focus on the non-relativistic chemical potential p = p — Mbaryon because the 
relevant scale of p would be much smaller than the baryon mass.) 

Thus, through the remainder of this section, we are going to impose a compres¬ 
sion force F on the multi-instanton system, assuming that the instantons form a 
uniform zigzag of some lattice spacing D and amplitude e, or a uniform straight 
chain for e = 0 , and vary D, e and the orientation moduli of the zigzag to seek 
the minimum of the chemical potential (118111 for any given combination of F and 
M3/M4. 

It is convenient to introduce a new variable £ = e/D to write the chemical 
potental as 


p = F x D + N C \MM% x D 2 p 


n 2 N c F m (C M 3 /M 4 ) 
20AM X D 2 


(182) 


where we have plugged the value of the function F, computed in appendix A of 
Ref. [30], minimized with respect to the orientation angles <p r , 0 


Pm,(Ci M 3 /M 4 ) = min^ r (0„; (£ = e/D), M 3 /M 4 ), (183) 


To find the local extrema of p with respect to £ and D one has to solve the equations 


aF m \ (dF m \ 

dDJz ’ 1 ) D 


0 . 


(184) 


As expected Eqs. (|184|) have two branches of solutions, corresponding to the 
straight chain e = 0 and the zigzag e > 0. However, for each branch, one may have 


h In fact there is no explicit formula for the M 3 /M 4 ) as the minimization of T with respect 

to the cj) angle has to be done numerically. Nevertheles s, i t is fairly easy to calculate numerically 
both the J~m function itself and its derivative dTm/d^W^ 
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Fig. 7. Phase diagrams of the zigzag and straight-chain phases in the compression vs. M 3 /M 4 
plane, F in units of N c yj AMM|). The straight-chain phase is colored pink, while four other 
colors — red, yellow, blue, and green — denote zigzag phases with different instanton orientation 
patterns. The first-order transition between phases are indicated by black lines, the second-order 
transitions by white lines. 


more than one solution for the lattice moduli (£,.D). In such a case global minima 
must be selected upon the evaluation of ft for each solution. Physically the presence 
of multiple solutions signals that the system is close to a first order phase transition. 
Indeed, varying the compression force F in this region one would typically find that 
the global minimum jumps from one lattice geometry to another. 

Altogether, there is seven different transition sequences for different M 3 /M 4 
ratios. The transitions can be detected from figures [7] and HI which depict the phase 
diagram of all the zigzag and straight-chain phases in two different planes!-^ the 
compression force F versus M 3 /M 4 and the linear instanton density p = 1/D versus 
M 3 /M 4 . To complete this section let us summarize a few particularly noteworthy 
features of these diagrams: 

• Since we assume M 2 <C M 3 , the straight-chain phase always has the anti¬ 
ferromagnetic order of the instantons’ orientations. 

• The very first transition from the straight chain to a zigzag could be either first- 
order or second-order, depending on the M 3 /M 4 ratio: for (M 3 /M 4 ) < 0.725 
the transition is first-order while for (M 3 /M 4 ) > 0.725 it’s second-order. This 
difference is due to different orientation phases of the zigzag immediately after 
the transition: for (M 3 /M 4 ) > 0.725 the zigzag has the same anti-ferromagnetic 
order as the straight chain, which allows a second-order transition; but for 
(M 3 /M 4 ) < 0.725 the zigzag has a different orientation pattern — the non- 
abelian NAl or NA2 — so the transition is first-order. 

• The non-abelian phases NAl and NA2 of the zigzag cover much larger areas of 
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Fig. 8 . Phase diagram of the zigzag and straight-chain phases in the linear density vs. M 3 /M 4 
plane; the density p = 1/D is in units of yAMMj. The stable straight-chain phase is colored pink, 
while the stable zigzag phases are colored red, yellow, blue, or green according to the instanton 
orientation pattern. Finally, the gray color denotes densities at which a uniform zigzag or straight 
chain would be mechanically or thermodynamically unstable. 


the phase diagrams on Figs. [7] and [5] than the abelian phases AF and AB. In par¬ 
ticular, at larger compression forces F — and hence larger chemical potentials 
/i, larger densities, and larger zigzag amplitudes — the instanton orientations 
usually prefer the non-abelian patterns. Only the backgrounds with M 3 ss M 4 
such as the model we have analyzed in Ref. [29] - favor the abelian orienta¬ 
tions. 

• Figure [ 8 ] has gray areas at which an instanton zigzag with a uniform lattice 
spacing D = 1/p and a uniform amplitude e (or a uniform straight chain for 
e = 0 ) would be unstable against instantons’ motion along the x 1 axis (the 
long direction of the zigzag). If we put A 1 instantons into a box of fixed 
length L = A/p and let them seek the lowest-energy configuration, they would 
organize themselves into domains of two different lattices with different lattice 
spacings and different amplitudes. 

• The NA1 phase of the zigzag occupies two separate regions of the phase diagram 
separated by the region of the other non-abelian phase NA2. The phase tran¬ 
sition between the lower-left NA1 region and the NA2 region is second-order, 
while the transition between the upper-right NAl region and the NA2 is weakly 
first-order: the lattice spacing and the zigzag amplitude are discontinuous across 
the transition, but the discontinuity is very small and hard to see graphically 
on figure [H 

• Likewise, the transition between the upper-right region of the NAl phase and 
the abelian AB phases of the zigzag is weakly first-order. On the other hand, the 
transitions between the anti-ferromagnetic phases of the straight chain or zigzag 
and all the other zigzag phases is strongly first-order, with largish discontinuities 
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of the lattice spacing and even larger discontinuity of the zigzag amplitude. 

• However, for small M 3 /M 4 the discontinuity becomes small; for M 3 /M 4 = 0 
it vanishes altogether and the phase transition between the straight chain and 
the NA1 phase of the zigzag becomes second-order. This is OK because for 
M 3 /M 4 = 0 or rather for M 3 = M 2 <C M 4 — the instantons’ orientations in 
the straight-chain phase are no longer anti ferromagnetic but form the period- 
four Klein-group pattern similar to the NAl pattern of the zigzag. 


5.3. Beyond Zigzag 

Finally let us report some preliminary results of the study of crystalline structures 
beyond the zigzag. This section summarizes a research in progress, to be published 
soon.®3 

To study the two-dimensional instanton crystals beyond a single zigzag-shaped 
line, one can run a numeric simulation of a multi-instanton system, where the instan¬ 
tons’ locations and orientations are both allowed to evolve seeking the minimum of 
the net energy. To obtain a two-dimensional crystal with two large dimensions we set 
M3 ~ M 4 but M2 = 0, which allows the instantons to spread out in two dimensions 
(a,’i,a,’2) due to two-body forces (117511 and (11761) . Thanks to homogeneous scaling 
of these forces with distance the net energy scales like P 2 with two-dimensional 
density, so instead of varying the external pressure we have simply confined 1000+ 
instantons to a large fixed area of the simulator. 

A priori, we did not know what kind of lattice the instantons would form: a 
square lattice, a triangular lattice (tiled by equilateral triangles), or something less 
symmetric. To our surprise, the simulations produces different lattices for different 
M3/M4 ratios: 

• For high ratios 0.80 < (M3/M4) < 1 the instantons form a square lattice with 
anti-ferromagnetic orientations, as shown in figure l9l (page [56]) . 

• For low ratios 0 < (M3/M3) < 0.55 the instantons also form a square lattice, 
but with a very different orientation pattern shown in figure [10] (page 1561) . We 
call this orientation pattern non-abelian since unit translations of the lattice in 
different directions are accompanied by anticommuting SU (2) rotations of the 
instanton orientations. 

• For the medium ratios 0.60 < (M3/M4) < 0.75 the instantons form a triangular 
lattice (tiled by equilateral triangles), while the orientations form a non abelian 
pattern shown in figure |TT] (page 1571) . 

To be precise, the numeric simulations did not produce clean monocrystalline lat¬ 
tices such as shown on figures [9] [Tl] but rather polycrystals made of small randomly- 
oriented pieces of such lattices separated from each other by messy boundary lay¬ 
ers. Fortunately, the nature of the ideal monocrystalline lattice was quite clear 
from such polycrystals for most values of the (M 3 /M 4 ) ratio. However, near the 
first-order phase transitions between the square and the triangular lattices, i.e. for 
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y(») = ± 1 , 
y(@) = ± it 3 . 


Fig. 9. 2D anti-ferromagnetic square lattice of instantons. 



Fig. 10. 2D nonabelian square lattice of instantons. 


0.55 < (M 3 /M 4 ) < 0.60 and 0.75 < (M 3 /M 4 ) < 0.80, the simulations produce very 
confusing polycrystals with both square and triangular domains. Strictly speaking, 
they might have also contained domains of some other types we could not identify. 

To pinpoint the phase boundaries between the square and triangular lattices 
— and also to check for other lattice types — we have derived analytic formulae 
for energies (per instanton) of a rather large class of two-dimensional instanton 
crystals. Specifically, we allowed for simple lattices (one instanton per unit cell) of 
any geometry, where the orientations of the nearest neighbors are related by two 
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Fig. 11. 2D triangular lattice of instantons. 


fixed SXJ ( 2 ) twists, 

y(m+l,n) = y(m,n)xTi and, y(m,n + 1) = y(m, n)xT 2 , (185) 

with the same Tj,!^ € SU{ 2) for all lattice sites (m,n). This includes the abelian 
orientation patterns with T 1 T 2 = +T 2 T 1 and the non-abelian ones with T 1 T 2 = 
-T 2 Ti. 

Minimizing the net energy per instanton as a function of lattice geometry and 
the orientation twists Xi and T 2 , we found four distinct phases for different M 3 /M 4 
ratios: the three lattices we saw in the numerical simulations, plus the non-abelian 
rhombic phase which obtains in a narrow range of 0.768 < (M 3 /M 4 ) < 0.783. As 
shown on figure IT2l (page l58l) . this lattice has similar instanton orientations to the 
non-abelian square lattice, but the unit cell is deformed from a square to a rhombus 
with axis ratio about 2.6. The phase diagram for the four lattice types along the 
(M 3 /M 4 ) axis is shown in figure [T3] 

This completes our presentation of infinite two-dimensional crystals. The next 
subject concerns the thin instanton crystals , which are infinitely long in the x\ 
direction but have only a few layers in the x 2 . The transition between a one- and 
a two-dimensional crystals goes through a sequence of such thin crystals, and back 
in Ref. [30] and Sec. 15.21 above, we have assumed that the first step in this sequence 
away from a straight one-dimensional chain was the zigzag. Obviously, we need to 
check this assumption and to investigate the other thin crystal phases beyond the 
zigzag. 

The first tool for this task was again a numerical simulation, in which both 
instanton locations and their orientations were allowed to evolve as they would 
seek the net energy minimum. This time, we have turned on a small M 2 ^ 0 (but 
M 2 <C AX 3 . 4 ), so at low densities the instantons would form a one-dimensional chain 
along the X\ axis, but at higher densities they would spread out in the a ; 2 direction. 
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Fig. 13. Phases of 2D crystals for different anisotropy ratios M 3 /A/ 4 . 


We did not prejudice the manner of such spread-out but let the instantons find 
the best way. We allowed the instantons to move in both X 2 and X\ directions, 
but in the X\ we put a movable wall subject to a controlled pressure. Indeed, as 
discussed in Sec. 15.21 controlling the one-dimensional pressure F, rather than the 
one dimensional L , allows the instantons to form a uniform lattice for any F. 

A few thousands of such simulations for different pressures and (M 3 /M 4 ) ratios 
gave us a crude phase diagram of thin instanton crystals. Some features of this 
diagram were quite surprising: 

• The abelian zigzag phase does not exist! Its appearance near the top of the 
phase diagrams in Ref. f30l and in Figs. ED[ 8 ] in Sec. [5] is an artefact of presum¬ 
ing the instanton centers form a zigzag when they actually prefer a different 
configuration — a thin slice of a two-dimensional anti-ferromagnetic square 
lattice. 
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• Even beyond the zigzag, there are no abelian phases with twist angles <fi ^ n. 
The only abelian orientation pattern for an instanton crystal — thick or thin 
— is the antiferromagnet. 

• While there is a whole series of anti ferromagnetic thin crystals for (M 3 /M 4 ) > 
0.8, the anti ferromagnetic zigzag phase does not exist! Instead, there is anti 
ferromagnetic period-three wave phase 



At (M 3 /M 4 ) > 0.7, increasing pressure causes a second-order transition from 
a single straight line to this wave, followed by a first-order transition to two 
parallel lines of instantons with two-dimensional anti ferromagnetic order of 
orientations 


• On the other hand, the non-abclian zigzag phases NAl and NA2 do exist. 

Not so surprisingly, at pressures beyond the AF wave or NA zigzag phases, the 
instantons form n = 2 ,3,4,5,... parallel layers which act as a thin slice of an infinite 
two-dimensional lattice, with a similar pattern of orientations. For such slices of the 
AF square lattice, the square sides are always aligned with the long and short axes 
of the slice, although the squares themselves may be distorted into rectangles. On 
the other hand, for the NA square and triangular lattices, the squares or triangles 
do not look distorted, but the alignment of their sides is almost random, which 
makes for polycrystalline structures in the numeric simulations. 

To get a more precise phase diagram — and to make sure we are not overlooking 
some phases due to poor convergence — we have followed up the numeric simulations 
with analytical calculations of net energies (per instanton) for all the phases we 
have seen or guessed. We have allowed some tunable parameters, such as square or 
rhombic deformations for square lattices, or deformation and global SU(2) rotation 
of the orientation pattern for the triangular lattice. Optimizing such parameters for 
each combination of pressure and M 3 /M 4 ratio and always choosing the phase with 
the lowest enthalpy H = E + LF, and hence chemical potential fi = H/N (at zero 
temperature), we have obtained the phase diagrams shown on figures fTdlflfil loages 



To summarize, the preliminary analysis of the two-dimensional and quasi-one- 
dimensional lattices reported here demonstrates that the naive intuition about the 
preferred crystal structures at give densities does not always apply. Eventually one 
has to resort to numerical simulations to see which structures are actually taken by 
the instantons. Numerical simulations may also have convergence issues, especially 
close to the first-order phase transitions. Without appropriate care some crystalline 
phases may be overlooked. 
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M3/M4 



1 2 3 4 5 6 7 8 9 10 11 12 


Fig. 14. Phase diagram of ID —> 2D instanton crystals, chemical potential // (in units of N c M 2 ) 
versus M 3 /M 4 ratio. The phases are labeled as follows: 

1 line AF straight line 

wave AF period = 3 wave 

12,13, B4 AF square or rectangular lattice with 2,3,4 layers 

Z, Z' zigzag, nonabelian orientations NA1 and NA2 

□ 2, D4 NA square or slightly rectangular lattice with 2 or 4 layers 

<^3,^>4.<^5,<^6 NA tilted square or slightly rhombic lattice with 3, 4, 5, 6 layers 

A3, A4. A5 NA triangular lattice with 3, 4, 5 layers 

Q3,Q4 NA rhombic lattice ( axis ratio > 2), 3 or 4 layers 


6. Summary 

In this review we have summarized the existing results on the description of the cold 
dense phase of holographic baryons as instanton lattices in the generalized Sakai- 
Sugimoto model. The main feature of a finite density system of instantons in such a 
setup is the competition between the two kind of forces: the Coulomb repulsion of 
instantons and the restoring (gravity-like) force transverse to the natural instantons’ 
alignment. The latter force is the result of the existence of the transverse holographic 
dimension that appears in the holographic description. However, in the effective 
description that we used here this force was merely an external force akin to the 
one in the electromagnetic traps used in the cold atom experiments to create one- 
or two-dimensional arrays of atoms. 

The outcome of the competition of the repulsive and restoring forces is the ex- 
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M3/M4 



2.5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 22.5 25.0 


Fig. 15. Phase diagram of ID —> 2D instanton crystals, Pressure F along the x\ axis (in units of 
N c y/AMMf) versus M 3 /M 4 ratio. The phases are labeled similar to figure ITT 



Fig. 16. Phase diagram of ID —> 2D instanton crystals, Linear instanton density along the x\ 
axis (in units of \/XMM 2 ) versus M 3 /M 4 ratio. The phases are labeled similar to figure fl4l The 
gray areas correspond to non-uniform (gas+liquid-like) mixtures of phases. 


istence of different lattice structures depending on the density of instantons. The 
repulsion is minimized, when the instantons are far apart and the density is low. 
Meanwhile the instantons will occupy positions at the bottom of the potential well. 
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As was demonstrated in a simple example in section 13.41 when the density is in¬ 
creased the repulsion energy may become so high that they will prefer to climb 
up the well to minimize the repulsion. The dimensionality of the original lattice 
changes in such a case. Overall, increasing the density will eventually make the lat¬ 
tice thicker and thicker in the transverse direction of the restoring force. The change 
of lattice structures will typically occur in a first order phase transition, as most of 
the lattice structures will not transform smoothly into each other. 

The physical meaning of this kind of transitions should follow from the meaning 
of the transverse direction. Indeed, the holographic dimension has the interpretation 
of an energy scale. Growing size of the lattice in this dimension as the density 
is increa sed is resembling the emergence of the Fermi surface for finite density 
fermions.^ Note however, that there is no fermions in this picture. The baryons 
here are classical objects and the exclusion principle is realized by the hard Coulomb 
repulsion. 

In this work we considered one- and two-dimensional baryon lattices. Results 
for the two-dimensional lattices were presented in Sec. 15.31 They make part of a 
work in progress!-®^. The dependence of the lattice configuration on the anisotropy 
ratio parameter M 3 —> M 4 was presented in Fig. [13] Three-dimensional lattices 
were only briefly discussed in the point-charge approximation in Sec. 13.41 A more 
detailed investigation of the three-dimensional lattices will appear elsewhere. 

The details of the transitions between different phases depend very much on the 
details of the interaction forces. Those are quite non-trivial in the case of instantons. 
The leading order (in 1/A) solutions are the ADHM instantons parameterized by 
position, size and orientation in the SU ( N ) space. The force between them is the 
NLO effect and it depends non-trivially on the orientations. Moreover, many-body 
interactions contribute to the net interaction force. As a result, calculating and 
minimizing the energy of a multi-instanton configuration is a complicated problem. 

In section[4]we reviewed the known exact solutions for the instanton lattices. To 
be precise we were only interested in the expressions for the instanton density, or 
even a related quantity det L defined by Eqs. CD and d- The solutions include 
the straight instanton chain (1851) and the abelian zigzag (11071) . As was demonstrated 
by a later study, only the anti-ferromagnetic chain still have chances to remain the 
ground state, in the isotropic setup M 2 <C M 3 = M 4 (Sec. 15.21) . As the preliminary 
results of Ref. EH show, the abelian zigzag yields its place in the phase diagram to 
other phases ffigure 17771 in section lA3l) . 

The solutions corresponding to the simplest configurations, such as the straight 
chain and the abelian zigzag are already sufficiently complicated. The formulae for 
the instanton density are so bulky, that there is not much sense presenting them 
in the paper. More elaborate configurations require even more resources. Therefore 
it seems no other efficient way than working in some approximation. In section 15.11 
we have reminded the reader that in the small-instanton approximation a D, 
the forces between instantons can be reduced to just the two-body ones. This ap¬ 
proximation allows to implement the numerical simulations in the search for lowest 
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energy configurations. In Sec. 15.21 we have demonstrated how the phase diagram 
looks like if we assume that the zigzag is the first non-trivial configuration at higher 
density. However, the preliminary results of the studies in Ref.[63]reported in Sec. 15.31 
show that for certain values of the anisotropy ratio, the zigzag phase is skipped and 
the chain transforms directly into another configuration ("Figs IT4l 1161) . This example 
shows that our naive intuition about the phase structure does not always work. 

Rather than fixing a given pattern of the instanton positions, like a zigzag, or 
any other fixed geometry, a more accurate method would be to leave the positions 
as moduli and then minimize the net energy to find a preferred configuration. It 
is hopeless to achieve this analytically, but a numerical simulation, which allows 
instantons to select their equilibrium positions, would do. Such simulations have 
already been performed for (quasi-) one- and two-dimensional lattices. The results 
were reported in this review. The transition from the infinite three-dimensional to 
four-dimensional lattices is a subject of a future work. We hope that the results 
of that work will allow to compare the phase diagram with that of the skyrmion 
lattices. 
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